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abstract 


The  results  of  calculation-  of  the  performance  of  a  planar  p-i-n 
thenscphotovoitaic  ceil  are  reported.  A  computer  program  has  been 
written  which  automatically  optimizes  the  cell  geometry.  The  program 
allows  the  bulk  recombination  rate,  the  surface  recombination  velocity, 
tr.e  hole  and  electron  mobilities  and  the  intensity  and  spectrum  of 
radiation  illuminating  the  cell  to  be  chosen  arbitrarily.  The  perform¬ 
ance  of  germanium  p-i-n  cells  is  investigated  fbr  monochromatic  rad¬ 
iation  and  radiation  from  18" 3°K  erbium  oxide  and  black  body  sources. 

Experimental  work  leading  to  the  fabrication  of  interdigitated 
p-i-n  thermophotovoitaic  cells  ia  described.  The  results  of  tests 
evaluating  the  performance  cf  fabricated  cells  are  discussed.  A  series 
of  experiments wse  carried  out  to  determine  why  the  performance  of  these 
cells  was  below  theoretical  expectations. 
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INTRODUCTION 

This  is  a  fins!  report  on  the  work  which  was  performed  under  costoract 
number  DAAB07-70-C-0129.  This  report  contains  a  detailed  description  of 
the  work  performed  from  June  1,  1971  to  June  1,  1$J2  under  the  above  con¬ 
tract.  Previous  work  covering  the  period  from  June  1,  1970  to  June  1,  1971 
which  was  performed  under  the  same  contract  has  been  reported  in  the  annual 
report  TR  ECOM-Oi29-l.  This  report,  which  will  be  referred  to  as  Report 
No.  1^  was  published  in  August,  1971  with  the  control  symbol  OSD-1366. 

The  purposes  of  this  contract  were  to  theoretically  investigate 
germanium  p-i-n  planar  photovoltaic  c*lls,  to  determine  the  optimum 
design  of  such  cells  for  thermophotovoltaic  use  and  to  fabricate  german¬ 
ium  p-i-n  cells  which  represent  the  present  state  of  the  art.  Figure  1 
shows  a  cross  section  view  and  a  view  from  the  contact  side  of  a  planar 
p-i-n  structure. 

One  of  the  purposes  of  this  investigation  is  to  determine  if  the 
p-i-n  structure  has  a  significant  advantage  over  the  st  ndard  p-n  photovoltaic 
cell  in  thermophotovoltaic  applications.  In  order  to  do  this, the  perfor¬ 
mance  of  the  p-i-n  cell  had  to  be  analyzed  under  a  variety  of  operating 
conditions.  This  was  done  with  a  computer  program  which  would  allow  the 
computation  of  cell  performance  under  a  variety  of  specified  incident 
radiation  spectra  as  well  as  a  variety  of  materials  parameters  such  as 
surface  recombination  velocity,  bulk  recombination  rate  and  mobility. 

The  computer  program  was  described  in  Report  No.  I.1  During  the  last 
year  of  the  contract  this  program  has  been  expanded  to  allow  for  automatic 
optimization  of  the  device  dimensions  when  the  incident  radiation  spectrum, 
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surface  recombination  velocity  and  the  bulk  recombination  rate  are  specified. 
The  first  part  of  this  report  deals  with  the  results  of  this  optimization 
program. 

Three  distinct  types  of  radiation  spectra  were  selected  for  study  - 
monochromatic  radiation,  black  body  radiation  at  1873°K  and  radiation 
from  an  erbium  oxide  radiator  at  l873°K.  The  analysis  with  monochromatic 
radiation  is  of  little  practical  value,  but  it  provides  insight  into  the 
performance  of  the  cell  with  other  radiation  spectra.  The  analysis '.with 
the  erbium  oxide  radiator  is  of  interest,  because  the  spectral  emissivity 
of  erbium  oxide  is  relatively  large  for  photons  with  energies  near  the 
direct  gap  in  germanium  and  relatively  small  for  photons  with 
energies  well  above  or  well  below  the  direct  gap  energy.  Thus,  radiation 
from  an  erbium  oxide  source  can  be  converted  into  electrical  energy  more 
efficiently  than  radiation  from  a  black  body  source. 

The  second  part  of  this  report  deals  with  the  experimental  program. 

It  contains  a  description  of  the  procedure  used  to  fabricate  p-i-n 
photovoltaic  cells  and  of  the  investigations  that  were  carried  out  in 
connection  with  the  fabrication  proc  ;ss .  The  results  which  have  been 
obtained  to  date  on  the  performance  of  germanium  p-i-n  photovoltaic  cells 
are  discussed.  A  series  of  tests  Was  carried  out  to  determine  why  the 
performance  of  these  cell;  ic  below  theoretical  expectations. 

In  the  third  part  of  the  report  the  most  significant  theoretical 
and  experimental  results  are  briefly  discussed  and  recommendations  are 


made  for  future  work. 


I 
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I.  THEORETICAL  INVESTIGATIONS 

A.  Introduction 

In  Report  No.  1'^  the  results  of  the  program  to  optimize  the  design 
for  a  device  illunimated  by  black  body  radiation  were  described.  It  was 
clear  from  these  results  that  additional  optimization  calculations 
needed  to  be  performed  to  determine  the  effects  of  various  spectra  as  well 
as  the  effects  of  changes  in  surface  recombination  velocity  and  bulk 
recombination  rate  on  the  performance  of  the  device.  A  new  program  was 
written  which  utilized  essentially  the  same  approach  for  evaluating  the 

i 

performance  of  the  device  as  the  program  described  in  Report  No.  1  but  which 
also  automatically  optimized  the  dimensions  of  the  device.  This 
Automatic  Device  Optimization  Program  (ADOP)  is  described  in  the  appendix. 

B.  Lifetime  Considerations 

In  Report  No.  1  we  pointed  out  that  in  the  usual  operating  range 
for  a  p-i-n  photovoltaic  cell,  the  lifetime  of  the  excess  carriers  will 
be  determined  primarly  by  radiative  recombination  rather  than  by 
recombination  through  trap  sites.  In  Figure  2  we  show  a  curve  of  the 
excess  carrier  concentration  as  a  function  of  minority  carrier  life¬ 
time  for  which  the  rate  of  recombination  due  to  recombination  through 
trap  sites  is  equal  to  the  rate  of  recombination  due  to  radiative  recombin¬ 
ation.  The  curve  is  also  plotted  in  terms  of  open  circuit  voltage  which 
can  be  correlated  directly  to  excess  carrier  concentration.  If  the  open 
circuit  voltage  for  a  particular  minority  carrier  lifetime  lies  above 
the  curve,  the  dominant  recombination  mechanism  will  be  radiative 
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re combination.  If  the  open  circuit  voltage  lies  below  the  curve,  the 

dominant  recombination  mechanism  will  be  recombination  through  trap  sites. 

Since  the  minority  carrier  lifetime  for  the  starting  material  is  in  excess 
-3 

of  2  x  10  seconds,  it  appears  that  for  operating  conditions  which  would  normally 

*  IT  **3 

be  found  in  the  operation  of  a  photovoltaic  cell,  lo  electrons  cnr)  the  oacinent 
recombination  mechanism  would  be  radiative  recombination.  If,  however,  the 
minority  carrier  lifetime  in  the  bulk  is  significantly  reduced  during 
processing  (say  to  a  value  of  10  to  100  sec.),  then  the  dominant  recombi¬ 
nation  mechanism  would  not  be  radiative  recombination,  and  one  should  base 
lifetime  calculations  on  the  usual  Hali-Shockley-Reed  type  recombination 
statistics.  Since  the  purpose  of  using  intrinsic  material  is  to  obtain  a 
very  high  minority  carrier  lifetime  in  the  bulk,  we  have  assumed  for  all 
of  our  calculations  that  the  dc  iin'-at  recombination  mechanism  is  radiative 
.recomb inatio.  •  radiative  recombine tion  parameter  is  found  to  be  etiual 

to  2.8  x  10"“'’  c../  cec“A  from  fundamental  thermodynamic  arguments.  None  of 
the  devices,  v.’nich  are  described  in  chrpter  II,  were  operated  in  the  radia¬ 
tion  recombination  region  during  performance  tests,  performed  under  this 
r  jntract.  The  reasons  for  this  are  discussed  in  chapters  II  and  III. 


^  »  « .0  c~ .L  O  v' 

The  reoults  oo- ■  i.~u  .rom  the  ADQ1  program  are  summarized  i..  Tables  I,  I’., 
V  and  VI.  Th:  results  fox  -.r.ochr  critic  radiation  will  be  analyzed  in  con¬ 
siderable  detail.  iy  of  i.i»-  icJLusions  arawn  from  these  results  will  be 
applicable  to  these  obtained  for  b.V>ck  bodj  radiation  and  radiation 
from  an  erblur  o.vid  source. 

1.  iionochrcr-a*'  'ati  Ou 

The  r*  ■'form'  uco  «~f  the  w-i-n  coll  ht  o  been  analyzed  for  throe  different 
wavelengths  of  r  t  o^b  remit!  c  ra liatio.i.  The  wavelengths  were  chosen  so  that 
their  absorption  inrgthc  would  vary 


ov&r  a  "ange  of  values  which  would  be 


approximately  equal  to  realistic  values  of  device  thickness  (xq).  The 
three  wavelengths  chosen  were  1. 50  ^m,  1.54  um  and  1.60  um  with  respective 
absorption  lengths  of  2.3  um,  21  and  104  um.  An  absorbtion  length 
is  defined  as  the  distance  into  the  germanium  at  which  63$  of  the  incident 
radiation  has  been  absorbed. 

The  results  of  the  output  obtained  from  the  ADOP  program  for  mono¬ 
chromatic  radiation  are  summarized  in  Table  I,  and  r  Figure  3*  The  spacing 
between  contacts  (Li)  is  fixed  at  25  un*  The  maximum  current  collectable 
per  unit  area  of  illuminated  surface,  if  every  photon  capable  of  generating 
a  hole-electron  is  absorbed  (CBHDOPT),  is  fixed  at  5  amps/cm  .  The  following 
observations  can  be  made  from  the  results  in  Table  I: 

l)  When  the  excess  carrier  concentration  is  very  large,  the  optimum 
device  thickness  is  determined  primarily  by  two  conflicting  factors.  The 
device  should  be  as  thick  as  possible  in  order  to  absorb  as  much  of  the 
incident  radiation  as  possible.  On  the  other  hand,  as  the  device  gets 
thicker  the  volume  in  which  recombination  can  take  place  increases,  while 
the  relative  increase  in  the  number  of  photons  absorbed  tends  to  become 
smaller.  Therefore  the  excess  carrier  concentration  must  decrease.  This, 
of  course,  results  in  a  reduction  of  output  voltage. 

The  dependence  of  open  circuit  voltage  and  short  circuit  current 
on  device  thickness  is  shown  in  Table  II.  The  data  in  this  Table  is 
obtained  from  the  output  of  the  program  from  which  device  U  was  optimized. 
Note  that  as  X0  increases,  the  open  circuit  voltage  and  maximum  power 
point  voltages  decrease,  while  the  short  circuit  current  and  maximum  power 
point  current  increase. 

In  cases  where  VSB  and  VST  are  both  0  (devices  #1,  11  and  20),  the 
optimum  thickness  (X0)  is  roughly  four  absorption  lengths.  About  90$ 
of  the  incident  photons  are  absorbed  in  four  absorption  lengths. 
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TABUS  II 

Dependence  of  Output  Parameters  on  Device  Thickness 

X  K  i.54  y,m 


m  « 

50  ,Aa  LP 

=  75 

VSB  “  VST  “  0 

xo 

voc 

CUBDSC 

mXPOK 

CMAXPOW 

PMAXPGW 

w 

(volts) 

(amp  sef2) 

(volts) 

(amp  cm"2) 

(watt  cm-2) 

25.0 

-5186 

3.475 

.Wt-1 

3.274 

1.454 

37.5 

.5127 

4.158 

.4385 

3.913 

1.716 

50.0 

.5075 

4.535 

.4334 

4.263 

1.848 

62*  5 

.5029 

4.743 

.4289 

4.455 

1.911 

75.0 

.4988 

4.858 

,4249 

4.558 

1.937 

67*5 

.4951 

4.921 

.4213 

4.6l4 

1.944 

100 

.4916 

4.956 

.4181 

4.642 

1.941 

125 

.4861 

4.Q86 

.4124 

4.661 

1.922 

150 

.4814 

4.995 

.4077 

4.66l 

1.900 

175 

.4773 

4.997 

.4039 

4.656 

1.879 

200 

-4738 

4.996 

.3999 

4.648 

1.859 

The  FORTRAN  variables  given  in  this  table  hre  defined  in  Table  XIV 
The  data  is  obtained  free  the  output  of  the  program  from  which  device  11 
of  Table  I  was  optimized. 


Figure  3.  surface  recombination  velocity  -  ck  sec' 
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2)  In  cases  where  VST  and  VSB  are  both  0  (device  $1,  11  and  20), 
PMfvXPOW  increases  as  the  absorption  length  decreases.  This  is  because 
the  optimum  value  of  XO  is  smaller  and  therefore  the  volume  in  which 
excess  carriers  are  generated  is  smaller.  Hence,  the  msaber  of  carriers 
generated  per  unit  volume  is  greater,  and  the  open  circuit  voltage 
developed  at  the  contacts  is  greater.  This  conclusion  can  be  collaborated 
by  a  comparison  of  open  circuit  voltages  for  devices  1,  11  and  20.  Note 
that  the  relative  variation  in  open  circuit  voltage  is  much  greater  than 
the  relative  variation  in  short  circuit  current. 

A  second  advantage  resulting  from  larger  excess  carrier  concentrations 
is  a  reduction  in  the  resistivity  of  the  bulk  and,  therefore,  a  reduction 
in  ohmic  loss  due  to  the  flow  of  current  Between  contacts.  However,  this 
effect  is  almost  negligible  for  the  large  excess  carrier  concentrations 
in  devices  1,  11  and  20. 

3)  It  would  appear  that  the  optimum  contact  widths  IN  and  LP  exhibit 
a  strong  dependence  on  the  surface  recombination  velocity.  For  instance, 
in  the  case  of  devices  2  and  3>  a  small  increase  in  the  values  of  VSB 

and  VST  from  0  to  10  cm  sec^causes  tr.x*  optimum  value  of  LN  to  change 
from  31  (im  to  119  and  the  optimum  value  of  LP  to  change  from  51  jim  to 
201  ^m.  These  changes,  however,  are  not  significant;  for  even  though 
the  optimum  contact  sizes  change  greatly,  the  output  power  is  quite 
independent  of  contact  sizes  for  small  values  of  surface  recombination 
velocity.  This  is  shown  for  devices  2  aiid  3  in  Table  III.  These  devices 
are  illuminated  with  radiation  of  a  1. 50  u &  wavelength.  The  output  power 
shews  similar  lack  of  dependence  on  LN  aw'  V*  for  samples  illumxuated  with 
radiation  of  other  wavelengths  when  surf,  ce  recombination  velocities 


are  small. 


TABLE  HI 


■  Variation  in  Maximum  Power  Output  for  Different  Contact  Widths 
(Devices  2  and  3  from  Table  I) 

XO  =  50 


VBS 

VST 

LN 

LP 

PMA2P0W 

(cm  sec-i) 

(cm  sec”'1') 

(pB) 

(|A®) 

(watt  cm"2) 

0 

0 

* 

30.6 

* 

51.0 

2.049 

5 

100 

2.048 

25 

100 

2.048 

12.5 

50 

2.048 

25 

37-5 

2.049 

10 

10 

_  * 
118.7 

200.9* 

1.963 

150 

250 

1.962 

250 

1.960 

50 

200 

1.961 

150 

150 

1.962 

^Optimized  dimensions 


4)  As  VSB  increeses  while  VST  is  held  constant  at  0  (devices  9,  10, 

18,  19,  27  and  28),  there  is  a  tendency  for  LN  and  LP  to  increase.  The 
effect  of  increasing  IN  and  LP  while  holding  LI  constant  is  to  decrease 
the  fraction  of  the  unilluminated  surface  which  is  subject  to  surface 
recombination.  Thus,  fev<er  excess  carriers  are  lost  dug  to  recombination 
on  the  unilluminated  surface.  This  tendency  for  th;;  contact  size  to 
increase  is  counterbalanced  by  an  increase  in  tb j  series  resistance  in  the 
device  as  the  average  distance  between  the  center  of  the  p  contact  and  the 
center  of  the  n  contact  is  increased.  Thus,  as  the  recombination  velocity 
on  the  rear  surface  increases  from  a  value  of  zero  to  a  larger  value,  we 
observe  LN  and  LP  increasing  until  the  series  resistance  begins  to  dominate. 

5)  In  cases  where  VST  is  increases  while  VSB  is  held  at  0  (devices 
15-17  and  24-26)  there  is  a  tendency  for  the  device  thickness  to  increase 
slightly  until  VST  exceeds  1.00  cm/sec.  The  tendency  for  the  optimum  device 
thickness  to  increase  can  be  accounted  for  by  realizing  that  as  the  surface 
recombination  velocit;,  increases  the  relative  importance  of  bulk  recombina¬ 
tion  velocity  deceases.  As  XO  increases,  the  average  excess  carrier  densit 
decreases.  Thir.  in  tm*»>  decreases  the  total  recombination.  Increasing 

XO  also  incre-zses  the  number  of  carriers  generated,  although  the  effect  is 
quite  small  when  XO  is  greater  than  three  absorption  lengths.  As  the  excess 
carrier  concentration  decreases,  VOC  decreases.  Thus,  the  optimum  value 
of  XO  represents  the  best  compromise  between  these  competing  needs.  The 
optimum  value  of  XO  will  change  as  the  relative  importance  of  the  surface 
recombination  dianges. 

When  the  surface  recombination  velocity  exceeds  100  cm  sec-1,  the 
optimum  thickness  tends  to  decrease.  This  is  because  the  excess  carrier 
concentration  becomes  low  enough  to  invalidate  the  assumption  that  the  ' 
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excess  carrier  distribution  is  fairly  uniform  in  the  bulk  region.  Since 
most  of  the  excess  carriers  are  generated  vitaLa  >ne  absorption  length 
of  the  illuminated  surface  and  since  these  carriers  must  diffuse  to 
the  contacts  to  be  collected,  the  excess  carrier  concentration  in  the 
region  of  the  contacts  tends  to  be  significantly  less  than  the  average 
excess  carrier  concentration  in  the  bulk  region  when  va-re  is  an 
appreciable  current  flowing  through  the  contacts.  Thus  tue  maximum  power 
point  voltage  will  tend  to  be  significantly  reduced.  The  device  tends 
to  become  thinner  to  minimize  the  reduction  of  excess  carrier  concern. ret ion 
in  the  region  of  the  contacts  due  to  the  diffusion  gradient. 

6)  In  cases  where  VST  and  VSB  are  held  equal  (devices  2-5,  11-14 
and  20-23),  the  effects  of  surface  recombination  occur ing  ou  the  illuminated 
surface  only  and  on  the  unilluminated  surface  only  are  combiiied.  Thus, 

IK  and  LP  increase  as  VSB  increases  to  100  cm  secT^  and  (for  devices 
11-14  and  20-23)  XO  increases  as  VST  increases  to  100,  but  starts  to 
decrease  as  VST  increases  to  1000. 

2 •  R&di&tion  from  an  erbium  oxide  source 

The  results  obtained  from  the  AD0P  program  for  devices  illuminated 
with  radiation  from  an  l873°K  erbium  oxide  source  are  summarized  in  Tables 
IV  and  V.  The  current-voltage  curves  of  these  devices  are  plotted 

Figures  4endg.  in  general  the  observations  made  for  devices  illuminated 

) 

with  monochromatic  radiation  are  applicable  to  these  devices. 

For  an  idealized  device  (device  29)  the  conversion  efficiency  is 
14,7£.  If  VST  is  held  at  0,  very  little  degradation  of  output  power 
occurs  for  values  of  VSB  smaller  than  100  cm  sec 3  On  the  other  hand,  if 
VSB  is  held  at  0,  there  is  a  significant  degradation  in  the  output  power 
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as  VST  is  increased  from  0  to  100  cm  sec-1.  This  is  because  the  area  of 
the  illuminated  surface  is  much  greater  than  the  area  of  the  portion  of 
the  unilluminaued  surface  between  contacts. 

For  devices  where  the  surface  recombination  velocity  is  the  same  on 
both  the  illuminated  an  unilluminates  surfaces,  the  majority  of  the 
surface  recombinations  occur  on  the  illuminated  surface.  However,  it  is 
expected  that  during  processing  it  will  be  more  difficult  to  maintain  low 
surface  recombination  velocity  on  the  un illuminated  surface.  Therefore, 
the  total  number  of  recombinations  on  either  surface  should  be  of 
comparable  magnitude. 

Seme  of  the  results  of  Table  XV  have  been  plotted  in  Figure  5 .  As 
would  be  expected  the  output  power  increases  fester  than  a  linear  rate 
with  increasing  incident  radiation  intensity.  However,  because  the  radiativ 
recombination  increases  as  the  square  of  the  excess  carrier  concentration, 
the  rate  of  rise  in  output  power  with  increasing  intensity  is  not  as  rapid 
as  would  be  the  case  if  all  recombination  occurred  through  trap  sites. 

The  effect  of  recombination  through  trap  sites  is  shown  by  the  results 
summarized  in  Tabic  V  These  results  are  plotted  in  Figure  6.  It  is 
interesting  to  note  that  t.ie  open  circuit  voltages  of  c ,vices  49  and  50 
are  greatei  than  tlv  open  circuit  volt-ag1'  o**  device  29.  This  is  because  the 
recombination  rate  tends  to  vary  linearly  with  the  excess  carrier  concentra¬ 
tion  (m/T  )  rather  than  with  the  square  of  the  excess  carrier  concentration 
(rn  ).  Therefore,  us  the  device  becomes  thirncr,  the>relative  increase 
in  excess  carrier  concentration  and  open  circuit  voltage  will  be  greater. 
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Consequently,  the  optimum  thickness  will  tend  to  be  smaller  and  the  open 
circuit  voltage  will  tend  to  be  greater. 

3.  Radiation  from  a  black  body  source 

The  results  obtained  fran  the  ADOP  program  for  devices  illuminated 
with  radiation  from  an  l873°K  black  body  are  summarized  in  Table  VI.  In 
general,  for  the  same  value  of  CURBMAX,  the  optimum  device  dimensions 
occured  at  slightly  larger  values  for  black  body  radiation  than  for  erbium 
oxide  radiation.  The  conversion  efficiency  is  also  lower  because  a  larger 
portion  of  the  radiation  energy  is  in  unfavorable  portions  of  the  energy 
spectrum.  Device  52  (Table  VI)  has  a  conversion  efficiency  of  12.7$  when 
illuminated  with  l873°K  blackbody  radiation  while  device  29  (Table  IV) 

has  c  conversion  efficiency  of  14.7$  when  illuminated  with  1873$  Er2°3 

-2 

radiation.  In  both  cases  the  value  of  CURDOFT  is  5.0  amps  cm  . 
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II.  EXPERIMENTAL  INVESTIGATIONS 


A.  Introduction 

The  discussion  of  the  experimental  work  will  consist  of  the  following: 
(l)  a  brief  description  of  the  device  fab,  ic&tion  process;  (2)  a  detailed 
description  of  investigations  involving  the  fabrication  process;  (3)  a 
discussion  of  the  performance  of  devices  that  have  been  fabricated,  and 
(4)  a  listing  of  the  deficiencies  which  exist  in  devices  that  have  been 
fabricated  and  a  description  of  tests  that  have  been  carried  out  to 
isolate  these  deficiencies  and  measure  their  effects. 

B.  Description  of  Device  Fabrication  Process 

The  stepB  involved  in  fabricating  a  TFV  cell  are  given  below. 

1.  An  intrinsic  germanium  wafer  is  chemically  polished  on  both 
sides.2 

2.  i\ a  nsulatint,  film  is  deposited  on  both  sides  of  the  wafer. 

3.  A  pattern  for  the  formation  of  the  n  contact  is  etched  in  the 
insulator  on  one  side  of  the  wafer. 

4.  Ketal(s)  for  the  n  contact  is  vacuum  evaoprated  on  the  wafer. 

5.  The  metal(s)  is  etched  so  that  it  covers  only  the  portion  of  the 
wafer  on  which  the  insulator  has  been  etched.  (See  Figure  "f.) 

6.  A  pattern  for  the  formation  of  the  p  contact  is  etched  in  the 
insulator. 

7.  Metal(s)  for  the  p  contact  is  vacuum  evaporated  on  the  wafer. 
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8.  The  metal(s)  is  etched  so  that  it  covers  only  that  portion  of 
the  wafer  on  which  the  insulator  has  been  etched  for  the  p  contact. 

9.  The  contact  metals  are  alloyed  to  the  germanium.  (See  Figure  8.) 
During  the  course  of  the  experimental  work  some  devices  were 

fabricated  according  to  the  above  procedure.  A  number  of  variations  in 
this  procedure  were  investigated,  and  a  variety  of  insulating  films  and 
contact  metals  were  investigated. 
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C.  Experimental  Investigations  of  the 
Fabrication  Process  • 

1.  Insulator  deposition 

The  insulators  which  have  been  investigated  are  pyrolytically 
deposited  silicon  dioxide  and  vacuum  deposited  aluminum  oxide. 

Silicon  dioxide  deposition.--  It  is  desirable  to  carry  out  device 
processing  at  the  lowest  possible  temperature  in  order  to  prevent  diffusion 
of  lifetime  destroying  impurities  into  the  germanium  bulk.  Therefore, an 
investigation  was  carried  out  ir  order  to  determine  the  lowest  possible 
temperature  at  which  the  silicon  dioxide  could  be  deposited.  A  reasonably 
fast  deposition  (900  £/min)  cculd  be  obtained  at  a  temperature  as  low 
as  310°C. 

A  method  of  device  fabrication  was  proposed  in  Heport  No.  1^  in  which 
the  p  contact  was  formed  before  the  n  contact.  Aluminum  was  to  bo  used 
for  the  p  contact.  After  the  aluminum  was  deposited,  etched  and  alloyed, 
silicon  dioxide  vas  deposited  on  top  of  the  aluminum.  Windows  for  n 
contacts  were  then  etched,  in  the  silicon  dioxide,  and  metals  (gold  and 
antimony)  for  the  n  contact  were  su> -.equently  evaporated  and  alloyed. 

The  necessity  of  etching  the  n.  contact  metals  was  to  be  eliminated 
because  of  the  insulating  layer  between  the  alvininum  and  the  n  contact 
metals . 

The  results  of  attempts  to  fabricate  devices  with  the  above  described 

> 

procedure  were  ac  follows:  In  cases  where  the  silicon  dioxide  layer 
was  deposited  on  top  of  the  aluminum  contact  at  36o°C,  the  p  and  n  contacts 
were  invariably  shorted.  Furthermore, the  silicon  dioxide  invariably 
cracked  on  the  aluminum  cace  tho  thickness  of  the  silicon  dioxide  exceeded 
2000  A.  Devices  for  which  the  silicon  dioxide  was  deposited  at  425  C 


were  occasionally  workable. 
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Aluminum  oxide  deposition. —  An  investigation  was  carried  out  to 

determine  if  vacuum  evaporated  aluminum  oxide  could  be  used  as  an 

insulating  layer  between  contacts.  It  was  observed,  however,  that, 

while  aluminum  oxide  appeared  to  cover  the  top  of  the  aluminum  contact 

well,  it  did  not  cover  the  contact  edges. 

Photoresist  problems  with  SiO„  films.—  The  positive  photoresist 
* 

AZ-1350  is  generally  used  during  the  device  fabrication.  When  this 

photoresist  is  used  on  pyrolytically  deposited  silicon  dioxide,  it  is 

necessary  to  take  certain  precautions. 

First  of  all,  the  developer  which  is  normally  used  with  this  photoresist 

(AZ-1350  Developer  )  contains  an  abundance  of  sodixsn  ions.  Since  silicon 

dioxide  is  highly  permeable  to  sodium  ion%  and  since  it  is  feared  that 

these  ions  will  have  a  detrimental  effect  on  the  surface  recombination 

,  * 

velocity,  it  is  desirable  to  use  the  developer  XP -7110-1. 

Secondly,  AZ-1350  has  a  tendency  not  to  stick  to  pyrolytically 
deposited  silicon  dioxide,  and,  therefore, it  is  necessary  to  take  extra 
precautions  to  insure  adhesion. 

2.  Contact  formation 

In  order  to  avoid  excessive  leakage  of  minority  carriers  through  the 

’  \ 

contacts,  it  is  necessary  that  the  contacts  be  very  heavily  doped  and 
fairly  thick.  The  requirement  of  low  temperature  processing  rules  out 
the  use  of  high  temperature  diffusion  techniques  because  of  the  danger  of 
diffusion  of  lifetime  destroying  impurities  into  the  germanium  bulk. 

Heavily  doped  contacts  can  be  obtained  with  alloying  techniques,  but  in 
order  to  obtain  sufficiently  thick  regrowth  regions,  it  is  necessary 
either  to  alloy  at  high  temperatures  or  to  deposite  rather  thick  films 
# 

Mfg.  by  Shipley. 
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of'  contact’ metals  (l^mor  thicker).  Experimental  investigations  have 
shown  that  the  thickness  of  vacuum  deposited  films  is  limited  because 
of  the  effects  of  surface  tension  in  the  films.  The  increased  surface 
■  -  tension  of  thick  films  has 'two  adverse  effects:  (l)  There  is  a  tendency 
for  the  films  to  ball  up  or  develop  craters  during  the  alloying  process. 

(2)’ The  adhesion  of  the  films  is  poorer. 

n  contacts 

k 

Kim  has  investigat  ’  several  combinations  of  metals  for  forming 
alloyed  n  contacts  and  has  found  that  the  best  results  were  obtained 
when  gold  was  used  as  a  carrier  with  arsenic  and  antimony  being  used  as 
dopants.  For  this  reason  investigations  were  carried  out  in  which  these 
metals  were  to  be  used  to  form  alloyed  n  contacts.  It  was  discovered 
that  vacuum  deposited  films  of  gold  and  antimony  adhered  very  poorly 
to  germanium,  and,  therefore,  it  was  impossible  to  etch  these  films  with- 

•  j 

out  removing  them  from  the  germanium.  Because  of  this  adhesion  problem, 
several  attempts  were  made  to  build  a  device  without  having  to  etch  the 
gold-antimony  film.  The  first  of  these  attempts  is  described  on  page 

The  second  attempt  was  identical  to  the  first  attempt  except  that  the 
•position  of  silicon  dioxide  on  the  aluminum  contact  pattern  was  ommitted. 
Instead,  a  second  layer  of  aluminum  was  evaporated  and  etched  to  cover  the  first 
layer  of  aluminum.  After  the  gold-antimony  film  was  evaporated,  it  was 
alloyed  to  the  germanium  at  a  temperature  below  the  aluminum-germanium 
eutectic  .  It  was  hoped  that  the  unalloyed  aluminum  film  would  not  be 
affected  by  the  gold-antimony  layer  on  it  during  the  alloying  process. 

Then  the  alloyed  gold-antimony  film  (which  would  stick  to  the  germanium 
during  the  photo  etching  process)  could  be  etched  to  separate  the  n  contact 
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from  the  p  contact.  It  was  evident,  however,  that  the  gold-antimony  film 
had  interacted  with  aluminum  film  during  the  alloying  process. 

The  proposed  steps  in  the  third  attempt  to  fabricate  a  device  were 
as  follows i  (l)  Silicon  dioxide  would  be  deposited  on  both  sides  of  a 
germanium  wafer.  (2)  A  pattern  would  be  etched  in  the  silicon  dioxide 
on  one  side  of  the  wafer  through  which  the  n  contact  would  be  formed. 

(3)  Gold  and  antimony  would  be  evaporated  and  alloyed  to  form  the  n 
contact.  (4)  The  gold-antimony  film  would  be  etched  so  that  it  would 
remain  only  in  those  regions  of  the  device  on  which  the  n  contact  puttern 
had  been  etched  through  the  silicon  dioxide.  (5)  A  p  contact  pattern 
would  be  etched.  (6)  An  aluminum  film  would  be  deposited  and  etched  so 
that  it  would  cover  only  tV.e  p  contact  region.  (7}  The  aluminum  would 
be  alloyed  to  the  germanium. 

After  the  gold-antimony  film  had  been  alloyed  [step  (3)3,  however, 
large  areas  of  balled-up  goli-antimony  alloy  formed  on  the  silicon  dioxide. 
When  step  (4)  was  carried  out,  it  was  observed  that  the  photoresist  broke 
down  before  the  thick  balled-up  regions  could  be  completely  etched. 

An  unsuccessful  attempt  was  made  to  improve  adhesion  of  gold  to 
germanium  by  evaporating  a  mixture  of  gold  and  germanium  simultaneously 
on  a  germanium  wafer. 

Eventually  it  was  discovered  tnat  gold  would  adhere  extremely  well  to 

germanium  if  the  gold  were  evaporated  immediately  after  the  germanium 

* 

surface  had  been  subjected  to  an  argon  ion  etch. 

The  ion  etch  is  carried  out  at  an  argon  pressure  ranging  from  50 
microns  to  10  microns  while  the  vacuum  system  ^.s  pumping  down  the  pressure. 
The  duration  of  the  ion  etch  is  about  30  sec.  The  germanium  wafer  is  in 
good  electrical  contact  with  a  cathode  that  is  maintained  at  a  voltage 
of  600  volts  D.C. 
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Several  methods  of  "etching"  the  gold  antimony  films  were  investi¬ 
gated.  These  methods  are  described  below: 

Chemical  etching. —  While  chemical  etching  is  procedurally  complicated, 
it  invariably  produces  sharp  etches  and  can  be  used  with  relatively  thick 
vacuum  deposited  films. 

Ultrasonic  "etching1.1--  Although  gold  sticks  very  well  to  ion  etched 
germanium,  it  does  not  stick  at  all  to  silicon  dioxide.  Therefore  step 
5  of  the  procedure  outlined  in  section  B  (page  24 )  can  be  carried  out 
by  removing  the  gold-antimony  film  from  the  silicon  dioxide  with  a  rigorous 
ultrasonic  cleaning.  Although  this  method  oi  "etching"  is  procedurally 
simple,  it  has  two  disadvantages:  (l)  It  does  not  produce  as  sharp  a 
pattern  as  chemical  etching.  (2)  It  tends  to  become  less  satisfactory  as 
the  gold-antimony  film  becomes  thicker.  This  is  because  the  adhesion  of  the 
film  to  germanium  becomes  poorer. 

Photoresist  lift-off  "etching'.'--  A  third  method  of  "etching"  can  be 
carried  out  as  follows;  (l)  A  photoresist  pattern  is  developed  on  the 
wafer.  (2)  The  metal  film  is  deposited.  (3)  The  photoresist  and  the  metal 
film  on  top  of  it  are  strijpoi  from  the  wafer,  leaving  a  metal  pattern 
which  complements  the  photoresist  pattern. 

This  method  of  etching  has  three  disadvantages:  (l)  It  does  not 
produce  as  sharp  a  pattern  as  chemical  evening.  (2)  It  does  not  work 
for  films  that  are  too  thick.  (The  photoresist  won't  lift  off  thick  films.) 
(3)  It  is  inherently  dirtier  than  chemical  or  ultrasonic  etching. 

The  advantage  of  photoresist  lift-off  etching  is  that  the  problem 
of  etching  the  second  metal  contact  pattern  without  attacking  the  first 


is  eliminated. 
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A  device  was  built  using  lift-off  etching  techniques.  The  processing 
steps  were  as  follows:  (l)  Silicon  dioxide  was  deposited  on  both  sides 
of  a  wafer.  (2)  Both  p  and  n  contact  patterns  were  etched  in  the  wafer. 

(3)  A  photoresist  pattern  was  developed  which  left  only  the  p  contact 
region  uncovered,  (4)  P.n  aluminum  film  was  deposited.  (5)  The  photo¬ 
resist  and  the  aluminum  film  above  it  were  stripped  from  the  wafer.  (6)  A 
photoresist  pattern  was  developed  which  left  only  the  n  contact  region 
uncovered.  (7)  Gold  and  antimony  were  deposited  after  an  ion  etch.  (8)  The 
photoresist  and  the  gold-antimony  film  above  it  were  stripped  from  the 
wafer. 

Arsenic  doped  contacts.—  Investigations  were  cairied  out  in  which 
arsenic  was  used  with  gold  or  gold  and  antimony  to  form  alloyed  n  contacts. 
Two  noteworthy  observations  were  made:  (l)  Arsenic-gold  and  arsenic- 
antimony-gold  film  tended  to  stick  fair.ly  well  to  germanium,  provided  that 
the  arsemic  was  evaporated  first.  However,  the  films  did  not  stick  as 
well  as  gold  antimony  films  of  comparable  thickness  which  had  been  deposited 
after  an  argon  ion  etch.  (2)  Alloyed  arsenic-gold  films  tended  to  ball 
up  much  more  readily  than  arsenic-antimony-gold  films  of  comparable 
thickness. 

Other  combinations  of  metals  have  been  investigated  for  n  contacts. 
These  investigations  are  discussed  later. 

p  contacts 

Aluminum  has  been  thoroughly  investigated  as  material  fo.  forming 
p  contacts  because  of  its  very  high  solid  solid  solubility  in  germanium. 

Thick  films  of  aluminum  (of  the  order  of  10  have  been  investigated 
because  of  the  need  to  form  thick  germanium  regrowth  regions  during 
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the  alloying  process.  Because  of  the  long  time  required  to  etch  these 
films,  the  photoresist  would  break  down,  and  pin  holes  would  develop  in 
the  aluminum  film.  Experiments  were  conducted  with  various  etchants  and 
etching  techniques  to  eliminate  this  problem.  It  was  determined  that 
photoresists  would  hold  up  satisfactorily  if  the  film  was  etched  in 
warm  aluminum  etch*  (  approximately  40°C) . 

Experiments  have  beer,  conducted  to  determine  if  the  adherence  of 
aluminum  films  could  be  improved  by  an  ion  etch,  (although  aluminum 
films  generally  stick  well  without  ion  etches) .  The  results  indicate 
that  the  adherence  of  aluminum  i3  improved  slightly  by  an  ion  etch. 

Other  metals  and  combinations  of  metals  have  been  investigated  as 
materials  for  forming  p  contacts.  The  results  of  these  investigations 
as  well  as  the  results  of  more  detailed  investigations  of  aluminum  will 
be  discussed  in  the  following  section. 

Alloying 

An  extensive  series  of  investigations  was  carried  out  in  which 
various  metal  films  and  combinations  of  metal  films  were  tested  for 
possible  use  as  p  or  n  contacts.  The  metallic  films  were  alloyed  to 
the  germanium  at  various  tenperatures  and  were  angle-lapped  and  stained. 
The  quality  of  the  films  was  judged  on  the  basis  of  the  smoothness  of 
the  crystal  regrowth  region  in  the  germanium. 

Junction  staining  solution.—  Some  preliminary  tdsts  were  carried  out 
in  order  to  select  a  suitable  junction  staining  solution.  Aluminum  halls 
were  alloyed  to  germanium  substrates  to  form  thick  heavily  doped 
regrowth  regions  in  the  germanium.  The  substrates  were  lapped  through  the 
alloyed  region  and  stained  vith  various  solutions.  The  best  results  were 

The  composition  of  aluminum  etch  is  given  in  Table  VIII. 
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obtained  with  a  30  second  stain  in  a  solution  consisting  of  20  parts 
HP  to  one  part  HNOy 

Alloying  technique.—  All  metallic  films  were  vacuum  deposited  on 
chemically  polished  germanium  substrates  which  had  been  cleaned  ultra- 
sonically  in  methanol,  acetone  and  trichlorethylene  and  exposed  to  an 
argon  ion  etch. 

The  alloying  was  done  in  an  alloying  system  in  which  the  substrate 

was  placed  between  two  carbon  heater  strips.  The  heater  strips  could 

be  heated  independently,  and  thus  the  sample  could  be  alloyed  in  a  temperature 

gradient.  The  spacing  between  heater  strips  was  about  a  half  millimeter. 

The  sample  was  placed  on  the  lower  heater  strip  with  the  surface  on  which 

the  metallic  film  had  been  evaporated  feeing  upward.  The  temperature 

of  the  upper  strip  was  maintained  about  100°C  higher  than  that  of  the 
* 

lower  strip.  The  alloying  temperatures  given  in  the  following  discussions 
are  those  of  the  lower  substrates. 

The  samples  were  heated  at  a  rate  of  about  10°C  sec.1  A  cooling  rate 
of  about  5°  C*  sec1  was  begun  immediately  after  the  alloy  temperature  had 
been  reached.  All  alloying  was  done  in  hydrogen. 

Aluminum  films. —  Aluminum  films  approximately  one  micron  thick 
were  deposited  on  germanium  substrates.  The  films  were  alloyed  at 
temperatures  of  450°C,  5>00°0  and  700°C.  The  samples  were  angle-lapped 
and  stained.  In  the  samples  which  were  alloyed  at  >00°C  and  700°C,  the 
regrowth  regions  were  extremely  irregular,  and  there  was  no  perceptible 

According  to  EerNissey  the  regrovrth  region  in  the  germanium  substrate 
tends  to  be  smoother  when  the  temperature  of  the  alloy  melt  is  maintained 
above  that  of  the  germanium  substrate. 
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regrowth  in  the  sample  which  was  alloyed  at  450°C. 

There  was  also  a  tendency  for  craters  to  form  on  the  alloyed 
aluminum  film.  There  was  no  perceptable  re  growth  in  germanium  regions 
under  these  craters. 

Indium  films.--  An  indium  film  approximately  one  micron  thick  was 
deposited  on  a  germanium  substrate  and  alloyed  at  650°C.  The  film  balled 
up  badly  during  the  alloying  process,  and  regrowth  was  visible  only 
under  the  balled  up  regions. 

Aluminum  on  indium  films.—  Indium  films  (approximately  one  micron) 
followed  by  aluminum  films  (approximately  one  micron)  were  deposited  on 
germanium  substrates.  The  films  were  alloyed  at  650°C,  500°C,  450°C 
and  350°C.  The  re growth  region  of  the  film  alloyed  at  650°C  was  fairly 
smooth,  but  it  got  progressively  more  irregular  for  the  films  alloyed  at 
500°C  and  450°.  The  regrowth  region  was  quite  smooth  for  the  film 
alloyed  at  350°C.  This  was  probably  due  to  the  fact  that  the  alloying 
temperature  was  below  the  germanium-aluminum  eutectic  (424°C), and, there¬ 
fore#  the  aluminum  film  did  not  liquify  during  the  alloying  process. 

The  "cratering"  effect  described  in  the  discussion  of  aluminum  films 
also  o^cured  for  these  films.  The  craters  tended  to  get  progressively 
smaller  at  lower  temperatures  and  were  imperceptable  for  the  film  alloyed 
at  350°C. 

Antimony  on  gold  films.--  Gold  films  (approximately  one  micron) 
followed  by  antimony  films  (approximately  one  half  micron)  were  evaporated 
on  germanium  substrates.  The  films  were  alloyed  at  700°C  and  525°C. 

The  films  balled  up  badly  during  the  alloying  processes, and  regrowth  was 
visible  only  under  the  balled  u<  regions. 


Antimony  films.--  About  three  microns  of  antimony  were  evaporated 
on  a  germanium  substrate  and  alloyed  at  a  temperature  of  700°C.  The 
antimony  balled  up  very  badly  during  the  alloying  process,  and  no  regrcwth 
was  visible  in  the  germanium. 

Tin-antimony-tin  films.—  Some  samples  were  prepared  with  a  vacuum 
deposited  film  consisting  of  two  layers  of  tin  with  a  layer  of  antimony 
between  them.  The  first  and  second  layers  of  tin  were  approximately  (me 
half  and  one  micron  thick,  and  the  antimony  layer  was  approximately  one 
half  micron  thick.  The  samples  were  alloyed  at  k75°C  and  650°C.  The 
regrowth  regions  were  quite  smooth  and  regular. 

Conclusions. —  The  most  promising  junctions  for  p  contacts  were 
obtained  with  the  indiura-alunh num  films.  The  most  promising  junctions 
for  n  contacts  were  obtained  with  the  tin-antiraony-tin  films.  An  attempt 
was  made  to  fabricate  a  device  using  these  films.  It  was  discovered, 
however,  that  the  adhesion  of  the  films  was  not  good  enough. 

An  attempt  was  made  to  form  tin-antimony-tin  contacts  on  a  device 
by  alloying  the  film  before  etching  it.  Silicon  dioxide  was  deposited 
on  both  sides  of  a  wafer  and  was  etched  on  one  side  to  form  an  n  contact 
pattern.  A  tin-antiraony-tin  film  was  evaporated  and  alloyed.  But  during 
the  alloying  process,  the  film  balled  up  badly  on  the  silicon  dioxide,  and 
apparently  the  surface  tension  created  in  the  balled  up  regions  on  the 
silicon  dioxide  was  large  enough  to  cause  the  film  to  also  ball  up  in 
the  contact  region. 

Further  investigations  need  to  be  carried  out  to  determine  the 
maximum  thicknesses  of  indium-aluminum  fil'is  and  tin-antimony-tin  films 
that  will  stick  well  enough  to  grrmaniian  to  withstand  ultrasonic  cleaning 
and  photoprocess ing. 
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Contact  buildup 

Since  the  thermophotovoitaic  cell  is  designed  for  very  high  intensity 
illumination  ,  rather  large  currents  will  be  collected  at  the  contacts ; 
and,  consequently, it  is  necessary  that  the  resistance  in  the  contact 
fingers  be  quite  low.  Therefore,  investigations  were  carried  out  in  which 
the  resistance  of  the  contact  fingers  was  reduced  by  plating  relatively 
thick  (10-20  micron)  films  of  gold  on  the  contacts. 

Attempts  were  made  to  plate  gold  directly  on  the  alloyed  aluminum 
and  gold-antimony  contacts  of  a  completed  device.  The  gold  plating  on 
the  aluminum  did  not  stick  at  all.  The  gold  plating  on  the  gold-antimony 
contact  appeared  to  stick  reasonably  well,  but  the  plating  on  the  gold- 
antimony  fingers  bridged  over  to  the  aluminum  fingers  in  several  places. 

The  following  method  of  building  up  the  contacts  satisfactorily 
eliminated  the  bridging  over  problem  and  permitted  a  buildup  of  bcth 
aluminum  and  gold-antimony  contacts:  (l)  Very  thin  chrane  and  gold 
films  were  evaporated  over  the  contact  surface  of  the  device.  (2)  A 
relatively  thick  gold  film  was  electroplated  on  the  vacuum  deposited  gold 
film.  (3)  The  chrome-gold  film  was  then  etched  to  separate  the  P  and  N 
contacts.  The  chrome  film  was  evaporated  to  form  an  adhesive  layer  with 
the  aluminum;  and  the  gold  film  was  evaporated  subsequently  without 
breaking  vacuum  in  order  to  prevent  oxidation  of  the  chrome  film  and  thus 
insure  the  adhesion  of  the  electroplated  gold  film.  Jt  is  desirable  to 
build  up  the  contacts  with  electroplated  rather  than  evaporated  gold, 
because  thick  films  of  evaporated  gold  would  produce  an  appreciable 
amount  of  surface  tension.  This  surface  tension  would  cause  stresses 
in  the  germanium  which  could  significantly  reduce  lifetime. 


3«  Fabrication  problems 


Cgnpatability  of  etchants  and  processing  chemicals. —  A  variety  of 
chemicals  and  etchants  are  used  in  the  process  of  fabricating  a  device. 

Some  of  the  etchants  which  are  quite  satisfactory  for  etching  the  metal 
or  insulator  for  which  they  are  meant  will  attack  other  areas  of  the 
device.  The  fabrication  process  must  be  carried  out  in  such  a  way  that 
this  doesn’t  happen.  Table  VS  lists  some  of  these  chemicals  and  etchants 
and  their  effects  on  metals  and  insulators  that  are  used  in  the  fabrication 
process.  The  compositions  of  the  etchants  listed  in  Table  VI  are  given 
in  Table  Vlfi. 

Problems  with  gold  contacting  aluminum.-"  The  following  process  was 

used  to  fabricate  the  devices  which  were  evaluated  in  the  tests  in  section 

II-D:  (l)  Silicon  dioxide  was  deposited  on  both  sides  of  a  wafer.  (2)  A 

pattern  for  n  contacts  was  etched  in  the  silicon  dioxide.  (3)  A  gold- 
* 

antimony-gold  film  was  deposited  and  etched  to  form  the  n  contact.  (4)  A 
pattern  for  p  contacts  was  etched  in  the  silicon  dioxide.  (5)  An 
aluminum  film  was  deposited  and  etched  to  form  the  p  contact.  (6)  The 
contact  metals  were  alloyed  to  the  germanium. 

A  pi oblem  occured  in  step  ( 5) •  Normally  it  is  necessary  to  heat  the 
device  to  a  temperature  of  about  100°C  for  a  period  of  about  30  minutes 
during  the  photoetching  process.  This  heating  period  caused  the  aluminum 
to  react  with  the  gold.  A  method  of  photoprocessing  without  heating  the 
device  was  developed  in  order  to  eliminate  this  effect. 


The  film  consisted  of  two  layers  of  gold  with  a  layer  of  antimony 
between  them. 
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table  VE 


Effect  of  Chemicals  and  Etchants  on  Metals  and  Insulators 


Silicon  Dioxide  Etch 
Aluminum  Etch 
Gold  Etch 
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Antimony  Etch 
Chrome  Etch 
Nitric  Acid 
J-100 


Water 
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1 


4  4  2  1  2  1 

4  4  4  4  4  3 

2A 


1.  Etches  rapidly, 

1A.  Does  not  attack  aluminum  but  li f ts  of  the  aluminum  by 
attacking  the  germanium  surface  under  the  aluminum  film, 

2.  Etches  slowly 

2A.  Causes  flaking  of  film  after  several  days  immersion. 

3.  Has  a  discoloring  effect. 

4.  Has  no  apparent  effect. 


Chrome 
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TABLS  VIH 


Etchants  used  in  Device  Fabrication 


Silicon  Dioxide  Etch; 

104  mli  Hydrofluoric  acid 
454  ga  Ammonium  Fluoride 

682  ml  Water 

Aluminum  Etch: 

380  ml  Phosphoric  Acid 

15  ml  Nitric  Acid 

75  ml  Acetic  Acid 

25  ml  Water 

Gold  Etch: 

27  gm  Iodine 

60  ©a  Potasium  Iodj.de 
600  ml  Water 

Antimony  Etch: 

3  parts  Nitric  Acid 

1  parts  Hydrochloric  Acid 
6  parts  Water 

Chrome  Etch: 

27  grn  Cerium  Sulfate 

100  ml  Nitric  Acid 

400  ml  Water 
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4.  Heat  sink  fabrication  and  device  mounting 

A  photograph  of  a  finished  device  mounted  on  a  heat  sink  is  shown 
in  Figure  9.  The  interior  of  the  heat  sink  is  hollow.  Cooling  water 
enters  and  exits  through  the  openings  at  opposite  ends  of  the  heat  sink. 

The  portion  of  the  heat  sink  under  the  device  is  milled  to  a  thickness  of 
l/l6  inch.  The  milled  surface  is  left  unfinished  to  facilitate  the 
conduction  of  heat  between  the  surface  and  the  circulating  water.  The  heat 
sink  is  anodized  in  sulfuric  acid  to  form  an  insulating  layer  of  aluminum 
oxide  approximately  5  ^m  thick.  The  milled  out  portion  is  sealed  at  the 
bottom  with  an  aluminum  block  which  is  epoxied  on  with  stycast  2850FT* . 
Banana  jacks  are  soldered  on  two  copper  contact  pistes  which  are  epoxied  on 
opposite  sides  of  the  heat  sink  with  Stycast  285OFT*. 

The  device  is  epoxied  to  the  heat  sink  with  Stycast  2850FT*  used  with 
the  catylist  24  LV.  First  of  all,  the  device  is  mounted  (the  illuminated 
surface  facing  downward)  on  a  glass  slide  with  beeswax.  The  glass  slide 
provides  mechanical  support.  The  device  is  than  epoxied  to  the  heat  sink 
with  the  p  and  n  contact  pads  extended  over  opposite  edges  of  the  heat 
sink.  After  the  epoxy  has  cured  the  glass  slide  is  removed  by  melting 
the  beeswax,  and  the  beeswax  is  removed  with  trichlorethylene.  Finally 
the  contacts  on  the  device  are  soldered  to  the  contact  plates  with  Indalloy 
Solder  $8  used  with  Flux  x2. 

Manufactured  by  Eranerson  Cumming*  1 

-ft- 

Manufactured  by  the  Indium  Corporation  of  America. 
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Pin  Cell 


p-i-n  Cell  Mounted  on  Water 
Cooled  Heat  Sink 


5»  Evaluation  of  heat  sink  performance. 

The  following  experiment  war.  carried  out  to  estimate  the  degradation  in 

the  performance  of  heat-sinked  cells  due  to  the  heating  effect  of  the  9.95 
-2 

watt-cm  illumination  source:  A  shutter  was  inserted  between  the  cell  and 
the  illumination  source  to  prevent  heating  of  the  cell.  The  1-V  characteristics 
of  the  cell  were  monitored  on  a  curve  tracer.  A  series  of  traces  on  the  curve 
tracer  was  photographed  immediately  after  the  shutter  was  removed.  The 
effect  of  heating  could  be  determined  by  observing  the  degradation  in  the 
I-V  characteristics  on  the  photograph  of  the  traces. 

In  cases  where  there  was  no  fluid  circulating  through  the  heat  sink,  a 
very  significant  degradation  occured  in  the  cell  performance.  In  cases  where 
there  was  water  circulating  through  the  heat  sink,  however,  the  degradation 
was  practically  imrperceptable. 
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D.  Performance  of  Fabricated  Devices 

Some  devices  wore  fabricated  with  the  process  outlined  pages  24  and 
25.  The  open  circuit  voltages  of  each  device  was  measured  in  a  sample  holder 
in  which  the  illuminated  surface  was  faced  downward  on  a  glass  slide.  A 
tungsten  lamp  illuminated  the  sample  through  the  glass  slide  with  an 
incident  intensity  of  0.7  watts/cm2.  No  measurements  were  made  to  deter¬ 
mine  what  percentage  of  the  incident  light  was  absorbed  by  the  sample. 

The  results  of  the  open  circuit  voltage  measurements  are  shown  in 
Table  IX.  The  Table  also  lists  the  unique  characteristics  of  each 
device. 

Two  note  worthy  observations  can  be  made  from  the  data  in  Table  IX; 

(l)  The  quality  of  the  device  generally  improves  as  the  film  thickness 
(especially  the  thickness  of  the  n-contact  films)  is  increased.  (2)  The 
quality  of  devices  alloyed  at  temperatures  at  or  above  550°C  is  better  than 
that  of  those  alloyed  below  550°C. 

The  open  circuit  voltage  of  device  -#4  was  measured  again  with  the  0.7 

2  2 
watt/cm  tungsten  lamp  being  replaced  with  a  9.95  watt/cm  tungsten  lamp. 

The  measured  open  circuit  voltage  was  0.28  volts. 

Devices  #1,  4  and  6  in  Table  IX  were  mounted  on  heat  sinks. 

Because  of  problems  encountered  with  breakage  during  processing,  the  active 
area  of  device  il  was  reduced  to  0.6  cm  x  0.55  cm.  Problems  were  encountered 
with  the  adhesion  of  electroplated  gold  films  on  devices  $1  and  4.  As  a 
result  the  performance  of  these  devices  was  significantly  degraded  after 
they  were  mounted  on  the  heat  sinks. 

The  I-V  characteristics  of  device  *6  are  shown  in  HgurelO.  She 
2 

9.95  watt/cm  illumination  source  is  used.  The  short  circuit  current 
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Voltage  -  Volta 


Current  -  Amperes 
Figure  10,  Voltage  Curves  for 


Device  #6  (See  Table  IX) 


gure  11. 


Current  -  Voltage  Curve  for  Device  &£■ 
Reverse  Bias  Applied  (See  Tanio  XX). 
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obtained  is  240  rna.  Since  the  active  area  of  the  device  is  0.55  cm  x  0.30 

cm  (See  Table  IX  ),  the  short  circuit  current  collected  per  unit  area 
-2 

is  1.5  cm  . 

An  attempt  was  made  to  estimate  the  maximum  attainable  short  circuit 

current.  A  reverse  bias  was  applied  to  the  device,  while  it  was  illuminated 

by  the  9.95  watt  cm  source.  The  reverse  characteristics  appeared  to  approach 

an  asymptote#  (See  Figure  11.)  The  asymptote  intersected  the  current 

axis  at  550  ma.  This  current  is  presumably  the  maximum  collectable 

current  of  the  device.  The  current  collectable  per  unit  area  for  the  0.55 

-2 

cm  x  0.30  cm  device  is  3«3  amp. cm  . 

E.  Investigation  of  Deficiencies  in  Devices 
The  deficiencies  in  the  performances  of  the  devices  that  have  been 
fabricated  can  be  attributed  to  any  combination  of  the  three  following 
factors:  (l)  excessive  bulk  recombination;  (2)  excessive  surface 
recombination;  (3)  faulty  contacts.  The  following  discussion  describes 
how  the  relative  importance  of  each  of  these  deficiencies  is  evaluated. 

1.  Contacts  deficiencies 

A  series  of  samples  was  prepared  according  to  the  following  procedure: 

(l)  Silicon  dioxide  was  deposited  on  both  sides  of  a  polished  intrinsic 
germanium  wafer.  (2)  A  pattern  of  dots  was  etched  through  the  silicon 
dioxide  on  one  side  of  the  wafer.  (3)  An  aluminum  film  was  deposited  over 


O 
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half  the  surface  on  which  the  dot  pattern  had  been  etched.: 

(4)  Metal  films  for  n  contacts  were  deposited  on  the  other  half  of  the 
wafer  on  which  the  dot  patterns  had  been  etched.  (5)  The  aluminum  and 
n  contact  films  were  etched  to  form  a  complementary  pattern  to  the  silicon 
dioxide  pattern.  (l.e. ,  the  metal  films  patterns  consisted  of  dots  which 
covered  the  areas  of  the  surface  on  which  the  silicon  dioxide  film  hod 
been  etched.)  (6)  The  aluminum  and  n  contact  films  were  alloyed 
simultaneously. 

The  open  cir  ,ii+  voltage  of  the  devices  fabricated  according  to  the 

above  procedure  was  ;.r.?asured  while  the  unprocessed  surface  of  the  device 

2 

was  illuminated  with  the  0.7  waT;t/cm  source.  The  testing  procedure  was 
exactly  the  same  as  that  used  to  obtain  the  results  shown  in  Table  IX. 

The  results  of  the  measurements  are  summarized  in  Table  X.  The  following 

conclusions  can  be  drawn  from  these  results:  (l)  In  general  the  open 
circuit  voltages obtained  for  devices  with  the  arsenic-doped  contacts 
is  better  than  that  obtained  for  devices  without  arsenic-doped  contacts. 
(2)  The  open  circuit  voltage  seems  to  hit  an  upper  limiting  value  of  0.20 
volts. 

2.  Lifetime  damage  due  to  silicon  dioxide  deposition 

A  series  of  tests  was  carried  out  to  evaluate  tn-*.  effect  of  silicon 
dioxide  deposition  on  the  bulk  and  surface  lifetime  of  devices.  Three 
germanium  strips  were  cut  from  an  intrinsic  wafer  which  had  been  polished 
on  both  sides.  Silicon  dioxide  was  deposited  on  one  aide  of  one  strip 
and  on  both  sides  of  another.  Four  gold  contacts  were  then  evaporated 
on  each  strip  to  form  the  type  of  structure  shown  in  Figure  12.  The 
portion  of  the  contact  on  the  rough  edge  was  for  the  purpose  of  forming 
a  good  ohmic  contact;  and  the  portion  of  the  contact  on  the  polished 
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TABLE  X 


Open  Circuit  Voltage  Output  of  Dot-Pattern  Devices 


Thickness 

Thickness 
of  n 

Device 

of 

Aluminum 

P^lm 

n 

Contact 

Contact 

Material 

Film 

00 

Alloying 

Open 

Circuit 

* 

.(A) 

Material* 

Temperature 

(°G) 

Voltage 

(volts) 

1 

4500 

Au 

360 

435 °c 

0.125 

Sb 

400 

Au 

120 

2 

4500 

Au 

360 

520°C 

0.17 

Sb 

*00 

Au 

XcO 

3 

4500 

As 

130 

435°C 

0.195 

Sb 

130 

Au 

700 

4 

4500 

As 

130 

520°C 

0.20 

Sb 

130 

Au 

700 

5 

4-500 

As 

200 

435°C 

0.20 

Sb 

500 

Au 

10,000 

6 

4500 

As 

200 

435°C 

0.15 

Gb 

500 

3 

Au 

3,000 

Sb 

20,000 

* 


The  a  contact  films  for 
were  evaporated „ 


a  device  are  listed  in  the  order  in  which  they 


Polished  Surface 


t  t  y  t 

Gold  Contacts  Rough 


Edge 

Figure  12.  Sample  Used  to  Test  Effect  of  Silicon 
Dioxide  Deposition  on  Lifetime 
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surface  was  for  making  contact  with  a  probe. 

The  sample  was  illuminated,  and  the  excess  carrier  concentration 
due  to  the  illumination  was  measured  by  measuring  the  conductivity  of 
the  substrate.  (See  Figure  13.)  Since  the  gold  contacts  were  not 
perfectly  ohmic,  the  conductivity  was  measured  by  a  four  point  probe 
technique.  A  fixed  current  flowed  through  the  outer  probes,  and  the 
voltage  drop  was  measured  across  the  inner  probes  with  a  high  impedance 
volt  meter.  Since  the  germanium  substrates  were  intrinsic,  the  ratio 
of  excess  carrier  concentration  to  intrinsic  carrier  concentration  was 
given  by  the  ratio  of  the  voltage  drop  without  illumination  to  the 
voltage  drop  with  illumination. 

The  results  of  the  measurements  are  summarized  in  Table  XT.  The  test 
fixture  and  illumination  source  are  the  same  as  those  used  to  obtain  the 
results  shown  in  Tables  IX  and  X. 

The  following  conclusions  can  be  drawn  from  the  results  in  Table  xi: 
(l)  The  deposition  of  silicon  dioxide  appears  to  degrade  the  bulk  and/ 
or  surface  lifetime.  (2)  The  silicon  dioxide  film  enhances  the  absorption 
of  useful  illumination.  (3)  The  effect  of  depositing  silicon  dioxide 
on  both  sides  of  the  wafer  is  not  significantly  different  from  the  effect 
of  depositing  it  on  one  side  only. 

3 •  Determination  of  the  dominant  source  of  recombination 

Theoretical  investigations  have  shown  that  surface  recombination 
velocity  is  smaller  for  accumulated  or  depleted  surfaces  than  for  flat 
band  surfaces.  An  experiment  was  carried  out  in  which  the  surface 
recombination  velocity  of  a  germanium  substrate  was  varied  by  varying 
the  surface  potential.  The  substrate  was  prepared  like  the  sample  shovn 
in  Figure  12  with  silicon  dioxide  films  deposited  on  both  sides. 


I 


Figure  13.  Apparatus  Used  to  Test  Effect  of  Silicon 
Dioxide  Deposition  on  Lifetime 


-51- 

TABLE  XL 


Effect  of  Silicon  Dioxide  Deposition  on  Excess  Carrier  Concentration 


Sample  4  Description 


Relative 
Excess  Carrier 
Concentration 
(n/nj 


1  IIo  silicon  dioxide  on  110 

either  side 

2  Silicon  dioxide  on  one  37 

side ;  side  without 

silicon  dioxide  was 
illuminated 

2  Silicon  dioxide  on  one  46 

side;  side  with 

silicon  dioxide  was 
illuminated 

3  Silicon  dioxide  on  both  40 

sides 
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Transparent  but  electrically  conducting  films  of  aluminum  were  deposited 
on  both  sides  of  the  substrate.  The  aluminum  films  were  electrically 
isolated  from  the  gold  contacts  and  «he  germanium  substrate. 

The  conductivity  of  the  samples  was  measured  with  the  apparatus 

2 

shown  in  Figure  13,  but  with  two  modifications:  (l)  The  0«7  watt  cm 
tungsten  source  was  replaced  with  a  low  intensity  source  in  order  to 
prevent  heating  the  substrate.  (2)  A  D.C.  bias  was  applied  to  the 
transparent  altaninuzn  plates  with  respect  to  one  of  the  voltage  contacts. 
The  results  of  the  measurements  are  shown  in  Table  XII. 

TABLE  XU 

Relative  Conductivity  as  a  Function  of  Surface  Bias 


Bias 

(volts) 

Relative 
excess  carrier 
concentration 

'  n.' 

3. 

10 

0.11 

5 

0.15 

0 

0.22 

-5 

0.24 

-15 

0.40 

-30 

0.49 

-41 

0.72 
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7'he  change  in  conductivity  resulting  from  the  surface  charge  layer  induced 
by  the  bias  was  appropriately  compensated  for.  The  following  observations 
can  be  made:  (l)  The  rate  of  recombination  on  the  surface  is  greater  than 
the  rate  of  recombination  in  the  bulk.  The  excess  carrier  concentration 
at  -41  volts  is  more  than  three  times  as  great  as  that  at  0  volts.  (2)  Con¬ 
trary  to  theoretical  predictions  the  excess  carrier  concentration  for 
accumulated  n-type  surfaces  (positive  bias)  does  not  increase  with  increasing 
bias.  Furthermore  the  excess  carrier  concentration  for  accumulated  p-type 
surfaces  (negative  bias)  does  not  increase  with  increasing  bias  as  fast  as 
it  Should.  Both  of  these  effects  are  probably  due  to  the  drift  of  ions  in 
the  unannealed  silicon  dioxide  film. 


4.  Conclusions 

The  following  conclusions  can  be  drawn  from  the  tests  which  were  carried 

out  to  investigate  the  deficiencies  of  the  devices,  (l)  For  devices  whose 

* 

open  circuit  output  voltage  is  well  below  0.20  volts  ,  the  dominant  deficiency 
is  faulty  contacts.  (2)  For  devices  whose  open  circuit  voltage  is  nearly 
equal  to  0.20  volts  ,  the  dominant  deficiency  is  excessive  surface  recombinations 
If  bulk  recombination  was  the  dominant  recombination  mechanism,  the  measure¬ 
ments  summarized  in  Table  Xllyould  have  been  independent  of  bias  voltage. 

These  conclusions  can  be  collaborated  with  a  mathematical  expression  which 
may  be  used  to  roughly  approximate  the  excess  carrier  concentration  at  the 
junctions  of  a  p-i-n  thermophotovoltaic  cell  in  terms  of  the  open  circuit  voltage: 


n  eV 

ni  "  exp  2kT 


For  V  =  0.20  volts  and  T  =  300°K,  n/n^  =>  47.  This  compares  reasonably  well 
with  the  values  of  46  and  4o  obtained  with  devices  2  and  3  in  Table  XJ. 


This  value  of  open  circuit  voltage  is  for  samples  illuminated  with  the 
0.7  watt  cm2  source  with  the  same  test  fixture  used  to  obtain  the  results 
in  Tables  IX,  X,  and  XI. 
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III.  CONCLUSIONS  AND  RECOMMENDATIONS 

As  ve  have  pointed  out  in  our  description  of  the  experimental 
results,  the  performance  of  the  devices  which  were  fabricated  falls 
far  short  of  the  performance  which  was  predicted  for  devices  with 
optimum  physical  parameters.  The  primary  reasons  for  this  are: 

1.  High  surface  recombination  velocity 

2.  High  bulk  recombination 

3.  Poor  junctions 

Any  future  program  which  is  directed  toward  the  improvement  of  p-i-n 
planar  photovoltaic  cells  should  take  these  points  into  consideration.  The 
surface  recombination  can  be  reduced  through  appropriate  annealing  procedures 
after  the  deposition  of  the  insulating  film.  The  alloying  problems  which 
have  been  discussed  can  be  solved  either  by  using  appropriate  alloying 
techniques  or  by  the  use  of  ion  implantation  to  form  the  n  and  p  contacts. 
The  problem  1th  bulk  recombination  .  eduction  during  processing  can  be 
eliminated  by  using  a  fabrication  procedure  which  keeps  the  maximum  process¬ 
ing  temperature  below  300 °C.  In  that  case  very  little  degradation  should 
occur  in  the  bulk  lifetime. 

A  maximum  conversion  efficiency  of  14.7$  was  computed  for  an  optimized 

« 

device  operating  with  illumination  from  an  l873°K  erbium  oxide  source. 

Nearly  half  the  energy  radiated  from  such  a  source  is  carried  by  photons 
whose  energy  is  below  that  of  the  germanium  band  gap.  The  use  of  a 
reflective  coating  which  will  return  these  photons  to  the  source  can  lead  to 
significant  improvement  in  the  operating  efficiency.  Future  work  should 
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include  the  fabrication  of  cells  which  have  such  reflective  coatings. 

Little  work  has  been  done  on  the  present  study  with  regard  to  heat 
sinking  of  the  device.  It  is  necessary  that  such  heat  sinking  be 
accomplished  because  of  the  very  rapid  degradation  of  the  device 
characteristics  with  an  increase  in  cell  temperature.  Heat  sinking  of 
a  p-i-n  planar  cell  is  more  difficult  than  for  a  standard  p-n  junction 
cell  because  of  the  fact  that  the  electrical  contacts  are  both  on  the 
seme  surface  as  the  heat  sink,  and  one  must  simultaneously  obtain  the 
electrical  isolation  and  good  thermal  contact  to  the  surface.  Any  future 
work  should  include  a  program  to  minimize  the  temperature  drop  across 
this  interface. 

If  the  above  mentioned  problems  can  be  solved,  and  all  indications 
are  that  they  can,  the  result  should  be  a  thermolphotovoltaic  cell  with 
relatively  hi$i  conversion  efficiency. 
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AFEENDIX 


In  this  appendix  the  ADOP  program  will  he  described  in  some  detail. 
The  input  parameters  and  the  printed  output  will  also  be  described. 

Description  of  the  Program 

The  program  is  described  very  briefly  in  the  flow  chart  shown  on 
page  58  •  First  of  all  the  input  values  are  set  and  preliminary 
computations  and  initializations  are  made  (Box  A) .  These  input  values 
include  an  initial  guess  of  the  dimensions  of  the  device  for  which  the 
maximum  power  will  occur.  The  input  current  is  initially  set  equal  to 
zero  (Box  B) .  After  additional  initial  computations  done  in  Boxes  B  and 
C,  the  open  circuit  voltage  is  evaluated  (Boxes  D,  E  and  F).  In  order  to 
evaluate  the  output  voltage,  a  boundary  value  problem  must  first  be 

■ft 

solved.  This  problem  is  solved  in  the  computer  program  by  successive 
over-relaxat 'on.  In  the  successive  ••'ver-relaxation  process  an  iteration 
(Box  D)  is  carried  out  in  which  a  new  and  better  estimate  of  the  solution 
of  the  boundary  value  problem  is  obtained  from  an  original  estimate.  Then 
a  convergence  test  is  carried  out  to  determine  if  the  new  estimate  is 
sufficiently  close  to  the  solution  of  the  boundary  value  problem.  If  not, 
the  new  estimate  is  used  in  the  iteration  process  to  obtain  a  still  better 
estimate.  By  repeated  iterations  and  convergence  tests,  the  estimate  is 
successively  improved  until  the  convergence  test  is  finally  passed. 

-  _________ 

The  theoretical  formulation  of  the  boundary  value  problem  and  a 
detailed  description  of  the  portion  of  the  computer  program  in  which  it 
is  solved  x.'  included  in  Report  Ho.  1. 
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adop  prosrm* 
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The  solution  of  the  boundary  value  problem  is  then  used  to  compute 
the  output  voltage  and  other  output  parameters  (Box  F). 

One*  the  open  circuit  voltage  is  computed,  a  new  value  of  current 
is  selected  (Box  H);and  the  computations  made  in  Boxes  C,  D,  E  and  F  are 
repeated  to  obtain  the  output  voltage  for  this  current.  This  process  of 
selecting  nev  values  of  current  and  evaluating  the  corresponding  values 
of  voltage  is  repeated  until  a  sufficient  number  of  values  are  obtained 
to  provide  very  good  estimates  of  the  short  circuit  current  and  the 
maximum  power. 

After  it  has  been  determined  that  the  short  circuit  current  and  the 
maximum  power  can  be  satisfactorily  estimated  (Box  G),  a  new  set  of 
dimensions  is  selected  for  the  device  (Box  J).  The  process  of  obtaining 
enough  points  on  the  I-V  curve  to  estimate  the  maximum  power  and  the  short 
circuit  current  is  repeated  for  the  new  set  of  dimensions. 

This  process  of  selecting  new  dimensions  and  evaluating  the  output 

parameters  for  these  dimensions  is  systematically  repeated.  The  dimensions 

are  selected  in  such  a  way  that  they  tend  to  he  close  to  those  dimensions 

•* 

for  which  maximum  power  will  cccur. 

After  it  has  been  determined  that  the  output  parameters  have  been 
evaluated  for  a  sufficient  number  of  sets  of  dimensions  (Bax  I) ,  the 
maximum  power  and  the  dimensions  for  which  maximum  power  will  occur 
are  determined  by  means  of  an  interpolation  process.  A  quadratic  function 
of  the  form 


*  ■  Ao  +  AA  +  A2*p  +  Vo  +  Vn  +  *?A  +  Vo  +  “w'n1* 

“13V0  + 


The  dimensions  which  are  varied  are  the  width  of  the  n  contact  (^n) » 
the  width  of  the  p  contact  {in),  and  the  thickness  of  the  device  (xc).  These 
dimensions  must  be  varied  in  discreet  increments.  (See  Report  No.  1.) 
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is  fitted  to  the  values  of  output  power  associated  with  ten  appropriately 
selected  sets  of  dimensions.  In  the  above  equation  P  is  the  power  out¬ 
put  and  the  A's  are  constants.  Note  that  since  there  are  ten  A's,  ten 
equations  (which  come  from  the  ten  sets  of  dimensions)  are  needed  to 
evaluate  these  A's.  Once  the  A's  are  evaluated,  the  values  of 
and  XQ  which  maximize  P  are  determined,  and  finally  the  maximum  power 
is  obtained. 

A  number  of  other  parameters  (such  as  open  circuit  voltage  and  short 
circuit  current)  are  obtained  from  quadratic  interpotations  for  the 
dimensions  at  which  maximum  power  occurs. 

Input  parameters 

The  input  parameters  are  listed  in  the  three  catogories  which  are 
defined  below:  physical  parameters— physical  values  which  ore  needed 
to  evaluate  the  performance  of  the  device;  dimension  parameters — numerical 
values  which  define  the  dimensions  of  the  device;  program  parameters — 
numerical  values  which  control  the  speed  and  accuracy  of  the  program  and 
the  amount  of  detail  in  the  output. 

Physical  parameters 

The  physical  parameters  which  are  listed  in  statements  1650-2100 
of  the  ADOP  program  are  described  in  Table  XIII  The  generation  rate  of 
hole-electron  pairs  as  a  function  of  distance  from  the  illuminated 
surface  is  read  in  at  statement  4250. 
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table  xiil 


LIST  OF 

PHYSICAL  PARAMETERS 

FORTRAN 

Physical 

Name 

Syabol(s) 

Comments 

XTOE 

kT/e 

(kT  over^e) 

DP 

D 

P 

hole  diffusion  coefficient 

DN 

Dn 

electron  diffusion  coefficient 

NI 

ni 

intrinsic  carrier  concentration 

RCOMFAC 

r 

recombination  factor 

VSB 

V6.B 

Surface  recombination  velocity  on 
‘bottom  or  unilluminataB  surface. 

VST 

V§T 

surface  recombination  velocity  on 
Top  or  illuminated  surface 

VGAP 

VG 

(energy  gap  voltage) 

CURDOPT 

(optimum  current  density) 

The  current  collected  pur  unit 
area  of  illuminated  surface  if 
every  photon  capable  of  generating 
e  hole-electron  pair  is  absorbed  and 
if  every  optically  generated  hole- 
electron  pair  is  collected. 


Dimension  parameters 

The  dimension  parameters  which  are  given  in  statements  3100-3400 
can  be  described  concisely  in  the  following  equations: 
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(1)  m  =  2  *  NLNl  *  K 

(2)  LP  =  2  *  NLP1  *  K 

(3)  XO  =  NXOl  *  H 

(^)  LI  =  NLI1  *  K 

(minimum  ^ 

(5)  4 value  off  =  (NXOKIH  +  l)  *  H 

(xo  j 

where 

LN  =  width  of  n  contact 

LP  =  width  of  p  contact 

XO  =  thickness  of  device 

The  fortran  variables  NLNl,  NLP1,  NXOl,  NLI1  and  NXOMIN  are  integers. 
The  variables  H  a«d  K  are  real  numbers  which  specify  increment  sizes  in 
centimeters.  (See  Report  No.  I.1)  The  dimensions  are  constrained  to  values 
which  are  integral  factors  of  these  increments.  In  the  ADOP  program,  the 
dimensions  of  the  device  are  varied  by  systematically  varying  the  values 
of  NLNl,  NLP1  and  NXOl  (Box  J  of  the  flow  chart) .  The  initial  guess  of 
the  dimensic  as  for  which  maximum  power  occurs  is  determined  by  the  initially 
specified  values  of  NLNl,  NLP1  and  NXOl. 

In  some  cases  the  optimum  thickness  of  a  device  will  be  much  thinner 
than  the  minimum  thickness  attainable  in  practice  due  to  limitations  in 
technology.  Therefore  the  program  has  been  adapted  to  optimize  the  device 
under  the  constraint  that  XO  will  have  a  specified  minimum  attainable 
value.  This  value  can  be  specified  by  appropriately  choosing  NXOMIN  and 
H  in  equation  (p). 

The  time  required  to  run  the  ADOP  program  can  be  reduced  by  choosing 
the  largest  possible  values  for  H  and  KL.  On  the  other  hand,  the  accuracy 
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of  results  obtained  increases  with  decreasing  values  of  H  and  K.1  Experi¬ 
mental  runs  were  carried  out  for  which  the  physical  parameters  and  the 
device  dimensions  were  held  constant  as  H  and  K  were  varied.  The  error 
in  the  value  of  maximum  power  was  made  negligible  by  choosing  sufficiently 
small  values  of  H  and  K.  The  error  in  the  value  of  maximum  power  obtained 
for.  larger  values  of  H  and  K  could  then  be  estimated  by  comparing  the 
maximum  power  obtained  with  larger  values  of  H  and  K  with  that  obtained 
with  the  small  values  of  H  and  X.  For  values  of  H  and  K  used  to  obtain 
the  results  in  Tables  I  through  VI  the  estimated  error  is  or  less. 

Program  parameters 

A  discussion  of  program  parameters  is  of  little  value  unless  one 
wishes  to  use  the  program.  In  this  case  a  detailed  description  of  the 
program  will  be  needed.  A  document  providing  such  a  description  will 
he  prepared  at  a  later  date. 

Printed  output 

Table  XIV  lists  and  describes  the  FORTRAN  variables  that  appear  in 
the  printed  output  in  the  order  in  which  they  appear. 

The  printout  can  be  divided  into  two  parts.  The  first  part  gives 
the  output  for  a  given  set  of  dimensions.  (See  Table  XV  )  This  part 
will  be  repeated  for  each  set  of  dimensions  for  which  the  points  on  an 
I-V  plot  are  to  be  obtained.  The  second  part  summarizes  the  outputs 
which  are  computed  for  each  set  of  dimensions  and  gives  information  which 
will  describe  the  output  of  the  optimized  device.  (See Table  XVI ) 

Output  for  a  given  set  of  dimensions.—  The  values  of  CURDNET, 
VOLTAGE  and  ETA  for  the  points  on  the  I-V  curve  are  listed  in  part  [B] 
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Fortran 

variable 

ILE 


PTKNT 

CURDNET 

VOLTAGE 

ETA 

voc 

VF 

CURDSC 

ETACOL 

PKAXPCW 


TABLE  XIV 

List  of  Output  Parameters 
Description 

(index  of  _1  inear  equations)  A  subroutine 
is  used  which  solves  a  simultaneous  system 
of  linear  equations.  If  ILE-0,  the 
system  is  solved  satisfactorily. 

Otherwise  the  system  is  singular  to  the 
accuracy  of  the  computer. 

(point  count)  an  index  which  is  chronolog 
ically  associated  with  a  point  on  the 
I-V  plot. 

(current  density  net)  net  current  out¬ 
put  per  unit  area  of  illuminated  surface, 
voltage  output  which  is  associated  with 
a  given  value  of  CURDNET. 
cell  effic^c.icy  associated  with  a  given 
pair  of  valves  of  VOLTAGE  and  CURDNET.  f.e 
ETA-  CURDNET  *  VOLTAGE/CURDMAX/VGAP. 

(open  circuit  voltage) 

(voltage  factor)  i.e,  VF  -  VOC/VGAP. 

(short  circuit  current  dendity)  The  short 
circuit  current  produced  per  unit  area  of 
illuminated  surface. 

(collection  efficiency)  i.e.  ETACOL 
-  CURDSC/CUREMAX. 

(maximum  power  point  jjower) 


0MAXPOW 


VMAXPOW 

CF 


ETACBLL 


XO 

LN 

LP 

LI 

K  J 

CURDOPT 


ISE 


(maximum  power  point  current)  T  he 
current  per  unit  crea  of  illuminated 
surface  at  the  maximum  power  point. 

(maximum  power  point  voltage) 

(curve  factor)  i.e.  CF  -PMAXPOW/VOC/ 
CURDSC. 

(coll  efficiency)  i.e.  EIACELL  -  PMAXPOW/ 
POWMAX. 

the  thickneta  of  the  photovoltaic  cell, 
width  of  the  n  contact 
width  of  the  p  contact 
spacing  between  contacts, 
increment  sizes 

optimum  cureent  density j  the  current  col* 
lected  per  unit  area  of  illuminated  surface 
if  every  photon  capable  of  generating  a 
hole-electron  pair  i*i  absorbed  and  if 
every  optically  generated  hole-electron- 
pair  is  collected. 

(index  of  symmetrical  equations)  A 
subroutine  is  uced^ which  solves  the  sim¬ 
ultaneous  system  of  equations  whose 
coefficient  matrix  mutt  be  positive 
definite  or  negative.  If  ISE*0,  the 
system  is  solved  satisfactory,.  Other¬ 
wise  the  coefficient  matrix  is  either 
singular  to  computer  accuracy  or  is  not 
positive  or  negative  definite. 
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'  i  i 


o 


ID 


me 

etace 

CURDMAX 


POWMAX 


ECAABS 


(indent  of  dimensions)  an  index  which  is 
chronologically  associated  with  a  set  of 
dimensions  for  which  the  points  for  an 
I-V  plot  are  obtained, 
same  as  EIACOL 
same  as  ETACELL 

(maximum  current  density)  the  current  out¬ 
put  per  unit  area  of  illuminated  cell  when 
all  generated  hole-electron  pairs  are 
collected. 

(maximum  attainable  povgr)  the  power 
output  per  unit  area  of  illuminated  surface 
when  all  the  generated  hole-electron  pairs 
are  collected  with  an  output  voltage  equal 
to  the  energy  gep  voltage,  I.e.  POHMAX" 
CUREMAX  *  VGAP. 

(CURDMAX/CURDOPT)  the  ratio  of  the  number 
of  carriers  generated  in  the  finite  think- 
ness  device  to  the  number  of  carriers  gen¬ 
erated  in  an  infinitly  thick  device. 
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TABLE  XV 


Printed  Output  for  a  Given  Set  of-  Dimensions 


ILE 


ILE 


[d3 

[E3 

[F] 

43 

1.318E+17 

1.318E+1/ 

44 

1.164E+17 

1.162E+17 

21 

9.904E+16 

9.873E+16 

20 

7.869E+16 

7 . 823E+16 

25 

5.249E+16 

5.187E+16 

34 

1.714E+16 

1.640E+16 

51 

5.758E+15 

4.981E+15 

=0 

45 

3.340E+16 

3.270B+16 

“0 

■=0 

54 

2.480E+15 

1.690E-5-15 

21 

1.379E+15 

5.823E+14 

30 

1.001E+15 

2.023E+14 

20 

8.708E+14 

7.079E+13 

20 

8.251E+14 

2.473E+13 

CA3 


S  : 


PTKNT 

CURDNET 

VOLTAGE 

ETA 

1 

0. 

4.439E-01 

0. 

2 

8.089E-01 

4.372E-01 

1.325E-01 

3 

1.618E+O0 

4.284E-01 

2.596E-01 

4 

2.427E+00 

4.159E-01 

3.780E-01 

5 

3.236E+00 

3.935E-01 

4.770E-01 

6 

3.892E+00 

3.286E-01 

4.791E-01 

7 

3.992E-K30 

2.548E-01 

3.8UE-01 

8 

3.654E+00 

3.680E-01 

5.038E-01 

n 

4.011E+00 

1.830E-01 

2.750E-01 

10 

4.016E+00 

1.171E-01 

* .762E-01 

11 

4.017E-K30 

6.494E-02 

9.774E-02 

12 

4.Q18E+QG 

2.676E-02 

4.028E-02 

13 

4.018E+00 

8. 595E-04 

1.294E-03 

VOC 

=  .4439 

VF 

=  .6726 

\ 

CURDSC 

=  4.018E-K)0 

'  1 

iJIACOL 

=  .9935 

1 

PMAXPOW 

=  1.345E+00 

| 

CMAXPOW 

=>  3.664E400 

1 

VMAXPOW 

=  .3670 

l 

CF 

=  .7539 

ETACELL 

=  .5038 

f 

XO 

=»  5.000E-02 

t 

LN 

=>  1.500E-02 

\ 

LP 

=>  2.500E-02 

\ 

LI 

=  2.500E-03 

1 

Cel 
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or  Table  XV.  The  column  labled  PTKNT  lists  indices  which  are  chronologi¬ 
cally  assigned  to  each  set  of  values  CLRDKET,  VOLTA®  and  ETA. 

In  part  [A]  the  reader  will  note  that  there  are  ten  rows  of  numbers 
excluding  the  rows  containing  the  expression,  ILE  =  0.  Each  of  these 
rows  corresponds  to  a  row  in  part  [B]  with  the  first  row  in  [A]  corresponding 
to, the  first  row  in  [B],  the  second  row  in  [A]  with  the  second  row' in 
[B],  etc.  The  cohan  designated  by  [D]  gives  the  number  of  iterations  required 
for  convergence.  The  edusn designated  by  [E]  gives  the  excess  carrier 
concentrations  in  the  middle  of  the  n  contact.  Thecdluan  designated  by  [F] 
gives  the  excess  carrier  concentration  in  the  middle  of  the  p  contact. 

The  expressions  ILE  =  0  result  from  a  curve  fitting  routine  which  is 
used  to  estimate  the  value  of  CURUNET  for  which  the  maximum  power  will 
occur.  (See  Table  XIV. . ) 

In  part  Cc]  the  output  parameters  and  dimensions  of  the  cell  are 
listed. 

Surmary  of  outputs  for  each  set  of  dimensions  and  output  of  optimized 
device.—  The  cell  dimensions  for  which  output  parameters  have  been 
computed  and  the  corresponding  output  parameters  are  summarized  in  part 
[D]  of  Table  XV  •  The  columns  headed  by  Ilf  and  LP  list  half  the  number  of 
increments  on  the  n  and  p  contacts,  and  the  column  headed  by  XO  lists 
the  number  of  increments  along  the  dimension  of  the  device  that  is  per¬ 
pendicular  to  the  illuminated  surface.  (The  values  listed  under  LH,  LP 
and  XO  are  equal  to  the  values  NLN1,  NLPJL  and  NX01  Ascribed  on  page 
62  ) •  Note  that  the  row  of  values  corresponding  to  the  number  9  in  the 
column  headed  by  II)  is  obtained  from  Table  XV. 

The  list  of  values  in  part  [ET  of  TableXVT  gives  the  dimensions  and 


output  parameters  of  the  optimally  designed  device. 
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The  array  of  points  shown  in  part  Cb]  of  Tableayi  may  be  represented 
algebraically  by  the  array  shown  below: 


X1 

X2 

X3 

vlin 

V1j£P 

vlxo 

V2in 

V2£p 

V2xo 

V3jfcn 

v3xo 

This  array  consists  of  the  eigenvalues  and  orthonoraal  eigenvectors  of  the 
symmetric  matrix  with  elements  A^.  The  power  may  be  expressed  in  the 
equation: 


^  .  vitn 
for  i°  1,  2  and 


+  \  v2  +  >  +  X 

A1  1  a2  v2  x3  3 

(n.M  -  n  )  +  v.  (n_ 
tP  iP'  IIP  jEP 


3. 


V 


+  v 


i*o(nxo 


Preceding  page  blank 


PROGRAM  MAIN  (  INPUT, OUTPUT,  TAPE5=IMPUT,  TAPE6=0UTPUT) 


7 


RE AL  'tN,  L"P,  L  I ,  ft  S' I N IT L  ,NO,X ,  M, MOP EO< ,  MAG IOUM  , MA  Gf  MA  X , N I", NOPPO',  N SMA X 
i  , KTOE »  NN,  tNP,  JJNCPOT , LNJOMJ,  NST2,  N02,  NS  4 ,NG4,KQ2DNErt, RU2DP6H  ,E  L  (200 
2i  »W1233),X(  200  )»Y(  200  1,  Z(  200),  CRN  T(  ICO),  VOLT!  ICO)  ,NSl  (22,49)  ,NS2  (2 
32  ,49  )  t  NS  3  {  22,  49  5,  MAGI  1(22*  49>,MAGI2  (22,  49) » MAG  13(22*4-9)  ,NS(22  ,49)  , 
4  MAGI (22,49), DUMMY { 22,  49),  DUMAG(  22, 49 ) , DNDT ( 22 , 49 ) »NEI  (22,49) ,CURDX 

5(49)  » CURDY  (49  j  _ _ 

“REAL  NS'8'ARVNSSBAR,"N'l2rG"BjLKX22T"  . 

REAL  NCINC,N1,N2»N3,NSRF 

INTEGER  R,S,  RH,  SH,  RPl,  $P  1,RHP  l,  SHP1 

INTEGER  RH2  . .  '  . "  ‘  " 

INTEGER  R I ,  R  I M I ,  S  IP  1,  SI 
J  NT  EGER  PATHKNT,  PATHNMB 

~TNrFGEirpmjr  ■  - - 

REAL  NLN,  NLP,  NXO,  NL  I,  VOCD(  50 ) , VFDt  50 ) ,  CURD  SCO  (  50)  ,E  TACOLD  ( 50)  , 

}  i’MXPOnl  D(  53  ) ,  CMX  POW  D(  50),VMXP3WD(  50) ,  CFO  (  50 )  ,E  TAG  (  50)  »CLRDMXD  (50) 

DIMENSION  XM<  3,1 ),  GIGVALt  3)  '  . . 

DIMENSION  CU  M  G  E  N  (  82) 

COMMON  IET,  JLNP,  JLNM,  JLPP,  JLPM,  JX3P,  J  XOM ,  LNLP  ,LNXO,  l.PXO  ,R  J  CLN  » 
i  RTCLPTRICX'O;  BM(3',  n;AMr3,‘3l7jMLNT50),‘jMLP(5C),JNX0(50)  ,Ii) 


LISTING  OF_PHY$ICAL  PARAMETERS 

I'.T  OE  =  D  .0259 
L’P  =  4  8.0 

[jN^96;'0'"  ~ 

.Jl  *2.5613 

pn=Ni  _  _ 

RC  OMF  A C=  2 . 8  E - 1 4  ‘  . . 

VS  B= LOO » 

VST»100.  _  _ 

\GAP=0 .65  . . 

CU  RDOPT  =  )  « 

l  STING  OF  PROGRAM  PARAMETERS" 


IC- 2 

><FSHNMB=1 
J  H=2 
LL=D 

I  MAX =150 
LCOPMI ft- 20 

iitvar  =1 

I  PAT H=  I . . 

KL!R\'P  =  0  5 

FPS I  LON- 0.001 

NC I NC=  5.0 

RF--0.35 

FACT  GR=  5.0 


200 
“  V'"- 
40y 
330 
5  33 
V  30 
300 
333 
3  53 
930 
1000 
1133 
1233 
1300 
5.303 
L3.!: 
13P-; 
J.350 
14: 3 
14  53 
15/0 
15'’  0 
IS'1  ) 
153:- 

r, 

J7  50 

1  a oo 

135  3 
1930 
1953 
200  3 

o 

210' 
21  "* 

2  2  00 
2  ?  5  2 
2303 
235-' 
?40‘i 
.244.1 

■)  ^  ">  'l 


J  ^  J 

25  5'' 

7  -t 


7  93  V' 


LISTING  OF  DIMENSION  PARAMETERS 
:  \r  ••  '3,  M*6 

■  r .  i_  _  t  j 

■ ci  -4 

’  o  i  'i~  2 
A  I ! -ft  l >-l 


■;  3 .  '■ 


i 


..  j  . 


LISTING  of  relaxation  factors 


OMEG  A8= l. 92  "  ~  “  . 

OF EGAS = 1.4 
Or'EGABH  =  1.92 

0FEGASB-0.9 . .  '  . . .  " 

0f£GAST=I.4 
OMEGAC® 1.4 
CFEGACH-1 a  4 

evaluation  OF  scale_factor  AND  READ  IS  AND  SCALING  OF  CUH3EN 

R£AD( 5  » 685  }  CJRDSCL 
S  C  ALE=  .CU  RDOP  T  /  C  J  RDS  CL 
DO' 3  NX =  I  f  B  2 
RE  AD  ( 5  » 685  I  CJMGEN'NX) 

CUMGENt NX)=CJMGEN(NX }*SCALE 
GSURF=  CUKGEN  (  2  )  ~CUHGEN(1 } 

LlSnNG_OF_  FIXED  CONSTANTS  AND  PRELIMINARY  COMPUTATIONS 
i  U=LMAX~l  H 

\LPHA=2.0*CN*DP/(  DN+DP) 

1  P=  Nl *N 1/ PC 
NOPPQ=  NC+PO 
OPNO=  PC+NO 
-C. ?N02=  POPNO+POPNO 
Li- PO+PO 
T,4=  P02  +  P02 
"r2=-NO<-NO 
iOt*  N02+ND2 

ECC= ( PC+NC5 / 2 ♦  _  _ 

\w  X  =  (  DP*PO+  L)N*Nti  )/  (  OF  «  QNT  . 

aE=N!+NI 
D=0 

t  r  =  i  c 

J3PT  S  =  0 
,DXX0M=3 _ 

•EGIN  LOOP  FOR  EVALUATION  OF  POINTS  ON  I-V  CURVE 

o  - :  d  t- 1 
LN=  NLN1 

an-  nlp  i 

v  0  '  NX  0 1  . ~  . ~ 

il  .',1  T  L=  NI.NI  +  NLP1+NL  I  i  » 

I  NI  T  L  =  NX  01 
;..N=  J,LN#K*2.D 
r  =  NLP^K^i .9 
r  fiL  I  ;< 

!.  ■  fJXO^H  -----  — . — 

uujNFr=a.o 
.  U=H/12.5E-4*0.5 

1  Of  r.  -  C‘J  MGF  N( \x  01  *N  I H+2  *  -  C  J’1  CEN  l  <  J 

,  ./a  - ; 

.  >'!' '  L  ,  A  -  J 

.'JHO,-*AX=  (  C,r,J  R  F  +  UJLKGEM  )  *  1 .6c-  19 
•OWMAX=Cj  RC^AX  >•  VC./P 
«' D :  <  LPt ;  N)/2  -  9<!  I 
•'  Ul)r.:  ALt-H/a  EPS  lLC)\/{  X3*X'0«V  OMO  ) 


3450 
3500 
3550 
35  0  0 
3650 
3700 
3750 
3900 
3950 
3900 
39  5  0 
4000 
4050 
4100 

41  50 

42  0  0 
...  5  0 
*333 
4350 
4400 
44  5  0 
4  500 
4  ..50 
4500 
46  6  0 
4700 
4 :  ->  :< 
<i  5  3  0 

'  It  *  '> 


45  k  ,i 
'•  \«  ,  0 
5  0  0  .. 
rvV; 
5i  0.J 

4  3  5.J 

5  200 


o '*:.o 


6r'4'> 


J. i  o  . 


J I 


5-;  y 

^45. 


cj’OOO 
605 '■ 

i  <\ 


f  ' 


1 ' 


•iOO  O  O  CJ  O'OOOO;  .  OOO 


i  V  i 

LISTING  OF  TERMS  THAT  MUST  BE  RESET  FDR  NEW  VALUES  OF  CURDNET  7  ’ 

_ _ _ _  7?* 

10"l=IINirL  “  “  7,/. 

J=  JI NI TL  .  '  ■  3u' 

MESHKNT  =  I _  _  _  31  : 

8  3. 

'  VALJE  OF  INOXITL  SET  T3  CAUSE  INITIALIZATION  35" 

a  V  ;* 

- lWiTL=5  '  "  '  '  ~  . . .  . .  39" 

9)  0 

EVALUATE  I  AND  J  DEPENDENT  TERMS  y?:, 

.  . ~ . .  "  . . . * '  '  '  . .  9}P 

COMPUTATION  OF  INCREMENT  SIZES  940 

95^ 

"STTREALTn - 

H= XO/ REAL  I  97  0 

REALJ=J  _  _  98  D 

~  K=  { (  LN+ LP) /2  • +L  I  )/REAL  J  ~  . .  . .  '  991) 

1000 

COMPUTATION  OF  INCREMENT  DEPENDENT  VARIABLES  1010 

.  . . . .  “  "  . . .  . .  ‘  102  0 

AK= ALPHA/ ( K*K  }  1030 

AH=  ALPHA/ (  H*H)  IDO 

”  H02DNE=H/2.0/DN/l  .6E-19 .  '  "  1050 

H02DPE=H/2.0/DP/1.6E-19  1050 

MDPFOK=-DP*l .6E-19/K  1070 

H2=~2.D‘>tt  *  .  *  "  “  *  "  . .  108  X 

DNE0K=DN*1.6E-19/K  10-0 

AH2PAK2  =  2.v,*AH+2.0*AK  1D> 

P.2  =2.0  <=RC0MFAC  1110.' 

HSOKS  =  H*H/K/K  U,v, 

HS  OKS  PI  =  HS  OKS  +1.0  HOO. 

HOATGS  =  H/ A}  PHA^GSURF  "  . .  . .  I  Wj  . 

DP2EK0H=DP*2.0n.60E-19*</H  117  V 

DN2EKOH=ON*?.U*1 «6E-19*K/ri  1130  ’ 

KC2DNEH=K/2.0/DN/ 1.6E-I9/H  '  11-Jl 

K02DPEH=K/2,0/OP/ 1.6E-19/H  l2->  ■  7 

H02 DNEK=  H/ 2 » 0 / DN/  1.6E-19/K  12  !?' 

.  HC2DPEK-H/2.0/DP/ K6E-19/<  . . . '  '  -  -  -  -  120'" 

HVSrO?A-H*V$r/2. 0/ALPHA  12  ,  -  - 

HVSBO?A=H*V$  6/2-0/ ALPHA  1 .20* 

C  121  ‘  « 

C  LISTING  OF  ITERATION  INDICES  12>' 

C  _  _  125- 

RH  =  LN/  (?V*K)  +0.5  ‘  "  “  .  *  l?Vi 

SHMLN/2.+-LI J/K+1.5  12" ! 

s=sh  i  ? :>  i 

RHPl  =  RH+  l 

R-  RHPi  l:> 

SH  PI  =  $  H  +  i  _____  .<  2  V  *>  ' 

'  "  RPi=R+l . '  "  '  .  '  '  ‘  '  "  1  -  JO. 

SP1-S  +  1  3  Hi 

JH-J+l 

jn  -  j<  2  r." 

IH= 1+2  ’  '  - 

1 N=  I  M  '  • 

JH Ml  =  Th-i  "  -  -  •  >  ,  ' 

..*W  =  JN-l  '' 

f  N ’-‘I  =  1  N- 1 

/HMJ  -  l  H~l  ' 


-7 


LISTING  OF  TERMS  DEPENDENT  ON  ITERATION  INDICES 

SPR=S-R 
RJNM25 JN-? 

rCRM=  Nl  f  (vs  r  -iVS8<:SMR/P,JNM2)/R2/X3 

TFRM2  =  V$T/R?/XU/R  JNM? 

TFRM3  =  VS  8/ R?  f  *  Q! RJNM2 
R'J  LKPT  S  -  (  IN-?  )*(  JN-2) 

EV  ALU  AT  I  UN  OF  G8JLK 

ICJ= I *NIH«2 
DO  Y)  NX=2,  INK! 

ICL= ICJ-NIH 

GG JLK( NX  I = ( C 3  VGEN (  I CJ ) -CJMGFN (  I  CL  )  ) /H 
ICJ= I  CL 

LIST  I  -1 G  OF  TFSTPUINT  INDICES 

INDXXC= (IN +2 )/2 
I NUXYC* ( JN+? i/2 
INDXXFs IN-1 
I  NUX  iNN=  II 
I  NDX  P.N=  S  PI 

LISTING  OF  TERMS  THAT  MUST  CHAN GE  WITH  FACH  NEW  VALUE  OF  ClKJNfl 

G  R  A  0  X  Q  R  =  {  G  S  J  R  F  ^  »  J  L  .<  G  F  N  -  C  J  R  0  N  E  T  /  1 . 6  T  -  1 9  )  /  X  0  /  -I  C  u  V  =  A  C 
NS  I  NIT  L---1  FRMHSwRI  l  I  rRMI*l  tHM  l+GBAJX.'H  ) 

INITIALIZATION  OF  MAGI  AND  NS 

IF  {  INOX  I!  L.NF.J)  GO  TD  '» t> 

INOX! T l- l 
IP?- 1+2 
JP?  =  J+? 

Of  2  3  UK- 1 , I p? 

DC  20  NY-  1 1  J  !’? 

*'A  G  [  (NX  » '  1 Y  J  i,v 
?  )  NS  (NX,  \n  }-  V>  !  NIT  i 

NS  '-'Ax-  F  mj  C  1  -  v  -  '<  i  * S  NR T  (  N  I  *N  1  ♦  GO. AN  >;K  )  ! 

GO JNUA AY  CCNlC,  I  A->  ON  MAGI  3N  T  l*»  A  A  C  N  SVr^.  !-  I;  -  *  . '» /» «.  i 

'i  ;  magi  “ *x -■  c  im  i  •  (  1 1  <i  v  }/  ?,  oml  i  ; 

JO  S3  ,  ;v 

5  J  MAGI  !  ?  »  S’  )  -- -:•!  /< 

OF  1  1  AT  l-.v  u‘^0  ‘  '  <•  >pr  « ?  +  3  ?*..<  >j  v'  ■* ft  j-  ■  ;  ) 

DFLTAl  C-U'-'l  GAC/i  -  i.  G  A  ‘  --L- rt 

HtSET  lCuPKNI  -  ; 

LTOPK NT • 3 


\  ■-  _ 
!  ■’  ' 
I  '.J 

3  6  I 
)  '!  ’ 
1  '* 

I  A 
!  '-3 
I  A  3 
!  ‘<4 
I  v. 
l'l '» 

1  •»  p 

U5 
;  S3 

1  A  j 
1  s  :> 

i  A  t 
I  A  *5 


1  s  ■. 


;  > 


i.  !  r  '  f  ['  r  '  ■ 


{  '  lG  ‘:r  .  :,f  .  ‘  !.  »  .  '.‘;f  .  '  ,  1 

,  ‘  ’Y  ;  N;  ,  _  0  J.'  . 

i  -'Mi  I  ,  s'  ">  J 
"■  7  )  N'  =  lil. 


>  '  I  » 


■  J  1 


-76- 

*0  WRITE  (6,550)  ( NS ' NX, NY  ) , NY=  1 »  JN  ) 

WRITE  (6,560) 

DO  8  0  NX  =  2 , T  H 

'  "  8  0  *  WR‘  f  T  E'  (6,550  )  1  MAG  I T  N  X ,~  N  Y  )“.  NY  =T »  J  H )  . .  ~  “ 

WRITE  (6,560) 

IF  {  UOOPKNT.  GE.LMAX  ).AND.(  INDXETA  .Ej3  .  1 )  )  GO  TO  228 
90  I F ( L0OPKNT .GE.LMAX  )  GO  TO’  515 
LOOPKNT=LCOPKWr+l 
C 

rr . P  "CENT "ACT --APPROX  IMAT  ION  'OF  MAGI  AND  NS,  AND  EVALUATION  OF  W  AND  V 

C 

IF  (SPI.EO.  JNM1)  NS ( 1 , S°  1  HNS  (  2,SPI) 
i+H02  0PEK*{MAGI(2,SHPl  )-MAGI(  2,SH)  ) 

IFfSPl.EQ.  JNM1)  GO  TO  115 
W  {  S  H )  =  0 . 3 

v  ( s  p  r )  =  i .  o  . *  .  . . 

DO  103  NY=SHP 1 , J  HM 1 
NS  T  2  =  NS  (2,  NY)  +NS  1  1,  NY) 

GP=DP2EK0H*( NST2  90PN0)/( NST2+N02) 

N$4=  NS  ( 2 ,  NY  + 1  )  +  NS  (  2,  NY  )  +NS(  1 ,  NY  + 1  )+NS(  1,NY) 

GN=  K02  DN£H*(  NS4+P04  )/ (  NS  4+P02PN0  2  ) 

'  PAGI  (2,NY1  =  maGI(2,NY-1)-GP>MmS(  2,  NY  )-NS(  1,  NYH 
VM  NY  )  =  W  (  NY-l  )+GPc-Y(NY  ) 

NS  ( 1 » NY  + 1 )  =  NS  l  l,.NY  )-6N*(  3  .0*(  MAG  I  (  3 ,  MY  )-MAG  I  ( 2  ,NY  ))  -  (  MAG  I  (*•  ,NY)--MA 
iGI  (3  ,NY)  )  )-(  NS  (2,  NY  +  i  J-NS  (  2, NY  )  ) 

100  Y  (  NY  + 1  )  =  Y  (  NY  )  *(  GN+GN+GN  )  *W(  NY  i 
C 

C  P  CONTACT— F/ALJAT  ION  OF  DELTA 

C 

NS  T  2  r  NS  (2,  JNMi  )  -.-MS  C  1,  JNM1) 

GP=DP2EK0H*l  NS  T  2  +POPNO  ) /(  NST2+N0  2) 

OELT  A=  (  PACK  2  ,  JII ) -MAG  l(  2,  JHMl  >+GP*(  NS(  2,  JNK1)-NS(  1 ,  JNH 1  )  >  /(W(..iMrM 
1 )+GP*Y( JNMl) i 


G  P  CONT  ACT— FV  ALJATION  OF  MAGI  AND  NS 

C 

DO  HO  NY  =  SHPl,  JHM1 
MAGI  ( ?  ,  NY  )  =  M  AG  I  (  2 ,  NY  )  H(NY  HDEtTA 
H  3  NS ( 1  , NY }  =  NS  (  l , NY ) n ( NY ) *DEL  TA 

NS  (1  f  JNPt  )  =  NS(  It  JNMI  )  +Y<  JNM  L  HOELTA 

r 

V 

C  NCCMT  ACT  —  APPRO<  I  MAT  ION  OF  MAGI  AND  NS,  AND  EVALUATION  OF  1,  A  NO 

C 

11:*  IF  (R.EU.2)  \a  il,  2  HNS  i  2,  2  J-.02DNEK*  (  MAG  I(  2,  2) -MAG  I  11,2)1 
IF  (  R «  E  Q » 2  )  GO  TO  135 
W  (  I  )  -  ) .  0 
Y  (  2 )  =  1  .  0 
00  123  NY  =  2 *  <H 
NST2-NS  (2,  NY  ) +NS  (  l,  NY  i 
6N- D N 2  f.  K tjH  •’(  N'5  T  2  *  PO ° N 0  *  /  (  N S T  2  •» P 0  2 ) 

NS  4  =  NS  (  ?.  i  \'Y  ♦  i  )  +  NS  (  2 »  ,Y ) +NS(  UNY+iHNSl  l, NY) 

GP-KU2 L>°r_h'  (  NS 4  »  isy.  )/<  NS4+^L2PNJ2) 


PAGI  ( 2  ,  NY  HM  AG  I  w  ,  NY  -  l  )  +  G\*<  N  S(  2,  NY  ) -N  S  (  1 ,  N  Y  )  ) 

WUiYHtf  (NY-!  )-GN*Y{  NY  ! 

NS  M  ,  NY  +  !  i  -  VS  (  I,  ,Y  )  *  l  i.C’M  MV*  II  3,  NY  J-MN  >  I  (  2  ,  N  Y )  )  -  (  r  AS  I  (A 

!  >  1  (  i  Y  )  )  >  -  i  v  >  (  /  *  . ,  i  :  N  s  (?,..}  i 
Y  (  ;.Y+  1  }  -  Y  ;  \  '  5  -  (  r  , 'r  > 


, !» N  ) 


Y 


1  /4 
I  H 
I  76 
i  7? 

)  73. 
17?. 

1  S3 

Ids 

167' 

1  \>  ?  ‘ 
13?' 
189. 
189' 
139: 

ho: 

191.. 

192  z 
m; 

j  9  4  .. 
195: 
196: 
19  73 
1  OD 
L  9  9  J 

20 i  , 
?o?;- 


2 » •  > . 
?(•-'• 
?  :>  3 : 
?d:o 
:  r  - 

m  r*  -» 

/ .  • 

::  *) 

F'i.O 
?n  '• 

;*  1 7 : 
2 1  *  :• 


•  t 


-  ,•  / 


N  CC"'-T;.rv--  t  At  J  A 1  I'.,N  Df  3LT  a 


NS  f?--  jS  (  7  »  *< !  -*  'o  (  1 ,  R  ) 


Ik 


:N=DN2C-K0H«{  NST  2  +POPNO ) / ( NST2+P3  2 ) 

OEir  A=  (HAGK  2,  RHP  1) -MAG  It  2,  RHJ-GNM  N S (  2  * ^  i-  NS(  1»R) ) )  /(W(RH)  -GN*Y(R 

1 )  )  _ _  .  _  ...  _ _  _  .  _  ... 

N  CONT ACT — EY  AL  JAT  IGN  Df  MAGI  AMD  MS 
DO  133  .NY=2 »  RH 

MAGI  12  ,  NY  )  =  M  A  G I  (  2 ,  NY  )  +W(  NY  )*DELTA 
133  NS  ( 1  »  NY  )  =  NS(  I »  NY  ) +Y (NY  )* DELTA 
NS(1,R)  =  NS(1,R)+Y{RI*DELTA 

NS  IN  INTERIOR  REGION 

15  DELT  AT  =  CELT  AT  C 
DO  150  NX  =  2, INM1 
"  DC  143  NY  =  2*  J  NM1 

DNDT  (NX,  NY  )=AH*(  NS!  MX  + 1 ,  \J  Y  )  -  N  S  (  N  X,  M  Y  ) -N  S  t  M  X  ,  N  Y  )  +  N  S  (  N  v-l  ,  NY)  }  ♦  4  K*  (  N 
IS ( NX , NY  +  1 ) -NS t  NX  ,  NY )-NS  (  N  X,  N Y ) +\ S( N X , N Y- 1 )) <GB  ULK { NX) -RC  OFF AC  *NS ( N 
2  X  ,  NY )  *  (  NS  (  NX  ,  NY  )  +NOPPO  ) 

NS  (NX  ,  NY)  =  NS  (  NX,  NY  )  +  ONDT  (  NX  ,  MY  )  *OEL  TA  T 
IF  ( NS ( NX, NY  ) .LT . 1.0)  N S(  NX , M Y )  =  1 .0 
140  IF  (NS (NX, NY ) .GE.NSMAX)  NS(  NX, NY )=N SVAX 
L-0  DF LT AT  =  DELT AI  B 

EV  ALJ  AT  ION  OF  NE  I 


DO  163  NX.  =  2,  INM1 

DO  163  NY  =  2 ,  J NM l  . 

1  6  J  NF  l  (  NX  .  \’Y  )  =  1 . 0  /  (  NS(  WX »  N  Y  5+ACCC) 

NAG!  IN  I  Ml  HR  I  On  REGION 


CM EG AH- CMEGACH 
DO  1  73  \X- 3,  I HMi 
MX  PI  *  NX  4  l 
NX  ML  a  MX- 1 
DO  1  73  NY  -  2,  JHMl 
MY  PL  NY  4-1 

MAGI  03  M=  (  (  NF  I  (  NX,  NY  PI  )+NE  I(  NX  ,  NY  M*MA3I  (NXP  l  -  N  V )  4-  (Nr.  I  l  NXM)  ,NYPJ  !  *  \ 
1  61  (NX Ml  ,  NY  )  )  4 ‘MG  l  (  NVM  J ,  NY  )  +  I  S  "K  $*  ( (  N £  !  <  OX,  N  Y*  1  )  +  NL  I  (  N  XH  ,NYH  )  I  A M A 
?G  I  (  NX  ,  NY  P 1)  +  (  N t  H  N  <  i  )  Y  )  +  N  (  I  (  N  XM  1  ♦  NY  )  I  '■  S‘AC  !  (  N  X ,  NY  •  1)  j  )  /  ( -(  SOI'  TRIM  ,F. 
31  (  NXfNYPl  )  <Nt"  I(  Nx«l,NYP  1  )  +V  j  (  'U,  -v  *.  '  C  <  >  i/  *l  r  NY'  )  » 

MAu !  I  ,'X  ,  NY  )  =  “  AG  1  I  NX  ,  NY  5  <  Jv  I  0/  w  0  \  u”  -  ~  v  A  -  1  »  N  f  ,  ’  v '  ' 

H  •  f.MtGAHs  (.Ml  vAiiH 

NS  AT  TCP  S3RF.CI 

DO  1B3  NY”  ?,  .1  Ntfl 

M-  NS  (  1  \,M  l ,  MY  )  +  H  )  A I  ON-  H  S I  3  ?A4  I  N  S  (  1  N  ■’  j  .  y  }  +  \>  5  (  |  , ,  \  Y )  I 
ISO  0$  U.N.NVUNS  I  IN,  NY  )  'G  •■'['GAS  1  •"  "—Nit  IN,  NY)) 


?i  .  i 
2  3.  ■  • 
2  3s  : 
2  3  1 3  ) 
24  3.)  N 
241 
2 42,  ii 
24  3  3  N 
24'.  3  ) 

246  >  ) 
2463  ) 

247  33 
2  4  S  3  J 
24933 
263  J  • 
2:  i 

7  s  J  •  ; 
2531  < 
2  5  V'  3 

n 

?■,!  1 
75 7  ;  > 
7  *5  -< .'  < 
>  6  ■; .  .. 
A  6.!  'Ml 
: 

r>  ■-  >  '• 


2  (?. 


ms  at  BfrriM  sjrfact 


i  r  s  rp  l  ,v»i  .n  5  v>  ■ .)  /ik. 

'•) .  /  1  ;  3  N  Y  "  ‘  I  •  ' 

M--  •  ........  (  N  M 

.  -  )  UN  (1  v:v  1  -  Nj  t  *  ,  •  ;  M  I  •'  5  i... 


.  i  - S  i  i .  '  i  ’  > 


NS  AT  IH"T  S  I  .A 


()'  Ml)  Nx-  :  *  )  N 


■  >  I 


M-  NS  (  NX  ,  2  )  “  . . . 

.NS  ( NX  , 1 J  =  NS  (  NX ,  1  )  +OM£GAS*(  M-N  S  ( N  X ,  i  )  ) 

MS  “A  r  ~R  I  GHr  ~E  dgF  " 

no  220  NX=1,  IN  _  _ 

M=  NS  I  NX  ,  JNMl  )  . . 

NS  (NX,  JN)  =  NS  (  NX,  JN)+0MEGAS*(  M-VS(NX,  JN  >  } 

' G V AT. U AT  I  ON”  OF  DEL T  ANS"  AND'  RESEt  OF  VALUE  DF'NS  . 

S  U  H=  0 «  0 

SUNS  =  0,0  ’  . . 

DO  222  NX~2  j  1 NH 1 
DC  222  NY = 2 , JNML 

5UM=$JM+NS(NX,NY  )  ""  *  '  '  . 

SU  MS  =  S  U  MG  +  NS  (  NX  ,  NY  )  *N$ ( NX , N Y  )  A/gM 

NS  B  A  R=  S  U  M/  fl  J  I.  K  P  r  S  /  y 

NSSbAR=SJMS/r  DJLXPFS  .  .  /£$/ 

S  t  GMANT=Q .0  A'  / 

DC  223  NY  =  2  ,  JNMl  / f^7 

S  r  GM'A NT  =  S  1  GM  ANf  +NS  (  IN* ,  NY  }  +NS*(  INM  1,  NY  )  *  '  . X°V/  . 

MGMANL=3»0  /A/ 

■F(RPl.GT.S)  GO  TO  225  Xy 

00  224  NY=.RPi,S 

? I GMANB=S ( GMANB+NS ( 2,  NY H-NSl  I, NY ) 

\G2  =  NS  3AR+T ERML 

Cs  NS  $  8 AR+N 12  *N$  BAR+T  ER  vi  2*  S  I  GM  AN  T  +  TER  M  3<“  SI&MANB-G3  AQXOR  . 

IF  1(002*802)  «LT.C  )  GO  TO  515 
)FLTANS  =  -  002  +  S ORT (  B02*B02  -  C  ) 

>0  226  NX  =  1,  IN 
?25  NY  =  l,  JN 
;S  ( NX ,  NY }  r  NS  (  NX ,  N Y  }  +0 EL T  AN S 

.  vSJRE  THAT  ENOjGh!  ITERATIONS^  HAVE  TAKEN  PLACE  FOR  STABILIZATION 
If  (LOOPKNT .LT .LOOPMINJ  GO  TO  60' 

0 MERGENCE  f  ES T 


?9?0 
203  0 
295  0 
3000 
301  0 
3  02  0 
303  0 
303  0 
303  0 

3031. 
303)  ( 

3032. 
303  2 1 
3032c 
303  3  - 
303  'A 
3034  1 
3  034* 
303«.  3 
303  .; 
303  O' 
33:;.  3 
3G3‘  • 
305-  ••  ^ 
303  /  ? 
3  0.3  i  * 
30330 
303-2 
A  J  ' .. 

3  0  ’  “  '* 
3  03  •? 

^  r  *>  i 
XOw  •  > 

3J..0 
3060  j 
30.'  <0 
i : 


■.•j  ODNS  =  A  PODS  *  NS  I  IN  I T  L 

F  (A»  UL’NS.IT  .  A8S(  CNCTI  2,  INDXNN  )  H  GO  TO  60 
l  F  (AEODNS  »U  .  APS  (  CNCTI  2,  IMDXPM  )  n  GO  TO  60 
‘  IA600\S-I  T  -  AOSl  f  NOT  I  INDXXF,  JNOXYC  I  )  )  OO  TO  60 
■  f  (Al-DONS.I  r  .  AIJSl  CNOT  (  INOXXG,  LNOXVC  ))  )  w0  TO  60 


3 ; . 

315 
?  i  0 

3  r-  ■ " 


i«R  I  EQJT  FOR  SJCCUSSFUl  EXIT 

(i  »rvAR.ro.n  co  m  25c 

•  Pi  1  T  t  (6,540) 

«  I  T  L  (&,5?0)  MESKNTtLOOPKNl 
1  .lire  i6,‘'40) 

!  0  2  30  .\/=  1  *  I« 

i‘‘  I  T  F  16,550)  :  Ni  (  \'.<,  NY  j,  NY-.  »,JN) 
•MTF  f6,,i.O! 
r  ?-. ,  \<  3,  i" 

*  »  r  T-  I  •!  M  N-. ,  NY  )  .  V;'r  i  ,  Jli) 

it '  ; 6 , 5 j 


k  }  )  }  *i 
\  7  }  n 

^  i  J  ^ 

12  30  3 

^  -X  f 

7  *\  “ 

325 


'■•JATIC”:  -I  TERMS  RilATFP  TO  3LM-DI  VOL  1  AGE  AMO  POWER 


_ _  -79- 

50  PATHNMB=  R-l  .  341 5'5 

34200 

..  Xp=SPl _ 1.34300 

KX=2  -  34400 

00  310  PATHKNT=1,  PATHNMB  34500 

_  _  _  _ _  _  _  34600 

EVALUATION  OF  EXCESS  CARRIER  CONCENTRATIONS  AT  JUNCTIONS  '  "  34700 

34800 

NN»(NS  12 , KN)  +NS  (  1,  KN )  )/ 2.0  34000 

NP=  f  NST  2 1 KP )  +NS(  1  *  k P  )  } /  2 * 6  “  . .  .  . . .  “  35000 

lFtIPATH.6C.il  GO  TO  255  35050 

WRITE  (6,580)  MES HKNT , P AT HKNT  _  ___  35100 

WRITE  (6,550  )  NN,  NP  '  “  .  . .  ""  35200 

’  CONTINUE  35230 

1  F  ( KN.  EQ.2 )  WR I T  E(  6,  610 )  L 30PKNT, NN , NP  35260 

~  fp  CTNN.GT  .070 )’.  ANCif  NPrGT.O.O)  )  "GO  T0”260  "  "  "  '  . . "'35300 

IF  (  (MESHKNT  .E0.M6SHNM8).AND.(<N.E3.2)  )  GO  TO  515  35400 

GO  TO  30  0  ________  _  35500 

. """  . . .  " . .  . .  -  —  -  -  . .  -  ■  35»£» 0 0 

INITIAL  POINTS  IN  THE  GRAPHICAL  INTEGRAL  ALONG  THE  X  AXIS  35700 

35B00 

>60  PTNE=  rN'N+NST27K"NV)/~2VO'+ACC(r  -  -  —  . . . 35500 

PT0HDY=(3.0*( MAG l ( 2, KN )-M AG  I  ( 2, KN- 1 ) ) /K  +  ( KaGI t 3 , KN ) -MAGI  (3,KN-i))/  36000 

IK1/4.0  36100 

GIXN=PT0HDY/PrNE*H/2.0  '  36200 

PT  NE=(NP+NS(  2,  KP)  J/2.0  +  ACCC  36300 

PT  DHOY~ ( 3oO*( MAGI ( 2,  KP ) — MAG  I ( 2,KP-1 1 )  /K  +  IMAG1 ( 3,KP)-MAGI  (3,KP-1 ))/  36400 

IK)/ 4  . 0  .  ’  . .  36500 

GIXP=PTDH0Y/PTNE*H/2.0  36630 

36730 

Evaluation"  of~graphical*  integraCs' "along*  th£' *x"  axi s . . 3&500 

36330 

IF  (KX.LT. 3)  GO  TO  280  37000 

DO  2 7 0”  NX  =  3,  KX  "  . .  "  "  ""  . .  37!  00 

PTNE= ( NS(NX-1,KN) +NS( NX, KN)  )/2.0+ACCC  37200 

PTDHDY=  (MAGI  (NX,  KN) -MAGI  (  NX,<N-l  U/K  37300 

Gl  XN=  GI XN+  PT  OHO/  /  PT NE*H  . .  . . .  “  . .  374  00 

PT  NE= I  NS ( NX-l *  KP ) +NS( NX, <P )  )/2.0+ACC:  37500 

PTDHDY= ( MAGI  (  MX , KP ) -M AG  I ( NX , <P- 1 ) )/K  37600 

27,j'*GtXP=GtXP+PT  DHDT/PTNE*H  '  37700 

37900 

EVALUATION  OF  GRAPHICAL  INTEGRAL  ALONG  THE  Y  A  X!  S  37900 

38000 

2CO  KP  Ml  =  KP-i  381  !>0 

GiY-0.0  382’>0 

*  0 0  2 9 0 ~  K Y  =  K N ,  K ? MY  "  "  '  . . .  38300 

PT NE=  (  NS  f  KX»  NY  )  +NS (  KX  »  NY  + 1  )  ) /  2.0+ACC-C  ,  38400 

P.TOHOX*  (MAGI  (  KX  +  l.NY  J-MAGK  KX,  NY  )  )/H  39600 

2 Gl  Y=GI  Y  +  PT  CHl)X/PTN6*K  33600 

38  7  3  3 

EVALUATION  OF  OJTPJT  VOLTAGE, _PDW6R  AND  EFFICIENCY  __  __  38300 

. . . '  .  38900 

OHMLOSS  =  KTOE/  l.6E-19.(  ON  +  DP  )*(  GtXN-GIY-  GIXP  )  *  39003 

DIFFPOT  =  KTOF*(  ON-DP  )/ t  DN  +  DP  )*ALOG(  (  N?+ACCC  )  /(NN+ACC:  )  )  39100 

JUNCPOT  =  KTCF*  At.  0G(  (  NN  «NO  )  *(  NP  *P3  1  /N  I  /N  I  )  3 '->2 00 

VOLT  AG 6=  J  J N C  & OT  +  01?  M>OT-3liM|.  OSS  39330 

PCWER=VCLI  AGE  *CJ  K  ON  FT  39407 

rr  a-  power/  poaImax  '  . "  “  39000 

8T ACJRD=CJRDMH  /CURCMAX  39630 

FTAV=VOLTAGE/VGAP  39700 

3  9300 
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WRIT  EOLit  OF  VaLJES  RELAT ED  T3  CURDNET,  VOLTAGE  AND  ‘POWER  "  3930 

•  4000 

IF  ( I  PATH. EC. 1 )  GO  TO  300  4005 

WRTTF16 , 5 50  )  CORONET',  VOLTAGE*  JUNCPQT,  dTFFPOT.'OHMLOSS  "  4010 

WRITE  (6,550  )  POW  6R,  PCMMAX,  ETA,  ETACLRD  ,  ETA  V  4020 

4030 

RESET  "OF* "PAT  H  OF  *  INTEGR  AT  ION  ~  '  "  "  . 4040 

4050 

0  KN=  KN-1  4050- 

TF  *  ( PAT'RNM  BCE  0C  n  *  GO  TOT  3T0  *  4070: 

KP  =  S  PI  +  (  JNM1  -SP1  )  *(  PATHXNT  )/(  PATHNMB- 1)  4030’ 

K:<  =  2+(  INM1-2  )  *(  PATHKNT)/(  PATHNMB-  L)  _  _  _  4030: 

0  CONT  I  NOE . . . . . .  41001 

4120r 

START  A  NEW  MESH  UNLESS  THE  FINEST  MESH  HAS  JUST  BEEN  COMPLETED  4130: 

'  ”  ~  41 4  DC 

IF  (KESHKNT.NE.ME5HNMB)  GO  TO  450  4130C 

41500 


BEGIN  STORAGE  AND  INTERPOLATION  WITH  PROCEEDURE'  FOR  FIRST  I  TERATION  41700 


IF  ( PTKNT.NF.  1)  GO  TO  340 

DO'  320'  NX  =  i'»~ !  N . 

00  320  NY  = 1 , J N 
o  NS2  .  NX  ,  NY )  =  NS ( NX ,  NY  ) 

00  330  NX  =  2 , IHM1 
DO  330  NY  =  2, J HM I 
•  MAGI2  (  NX, NY)  »MAGHNX,  NY  ) 

'  V2-  VOLTAGE - 

VOLT  1 1  )  =  V 2 
C2  -*C  JRONET 

I'.RIvT  (1  )  =  C2  . . . 

ii':'KNT  =  2 


41300 
41900 
42000 
42100 
42200 
42300 
42400 
42500 
42500 
4  27  00 
42300 
42900 
43000 


43020 

^VALUATION  OF  'CJRRENT  AND  VOLTAGE 'INCREMENTS  .  '  4  3  3  *>  j 

4  305? 

f.  1  NC  RE  G=  CJ  RDM  AX  /  N  C  INC  4  3030 

V  .’DINC-CINCRFG  . .  .  . .  43090 

CU RONE T  =  CJ RON f  T+CURD INC  4  3100 

GOTO  _10 _  _  _  _  .  _  43200 

> Nl ERPOLAf ION  PROCEEDUR6  FOR  SECOND  ITERATION  43400 

43503 

0  IF  (PTKNT*NE.2>  GO  TO  390  4  3500 

DO  ^50  NXri,  JN  43700 

OL  3  50  NY -  1 ,  J  N  4 33 00 

0  NSlfNX,  NY)=N$  (NX,  NY  )  ~  .  . . 43900 

DC  353  NX  =  2» IHMI  I  44003 

00  350  NY=  ?»  J HM1  44100 

"AGIl  (NX  i  NY)  sMAGK  NX,  NY  )  ”  '  44200 

VI  =  V()LT  AGE  44J00 

YULI  (2  i  =  VI  44400 

Cl  s'CJRDNET'  .  *'  '  -  '  '  *  4450  0 

ORIvT  (2  )  =  CI  44400 

PV KNT  =  3  44700 

cuRONF^cjon-Nn  *  cj^dinc  448?o 

G?  3  73  \r.s  1,  IN  .  449  0  I 

'JV  3  70  NV=UJN  44?0:: 

;.S  .'nX,NM4NS)(Nx,NYJ-N5?(  NX,'NY)  . . .  "45i  00 

jn  3,0  1  2f 1 HMl  ;52J0 

’■)(  '.3  0  !•.  4 'j33? 

>U  Gi  (  NX, NY )  =  MAGI  1 1 NX,  NY  )+MAGIi(NX«NY)-MAGI2(NX,NY)  4F4?0 
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GO  T 0~440 

BEGIN  QUADRATIC  INF  ERPOL  AT  IONS  WITH  IN  I TIAL I  ZA  TIONS 


39;  C3  =  C2 

C2  =  C1  _ 

Ci=CURDNEf 
CRNTIPTKNT  )=Cl 
V3  =  V2 
V2  =  V  1 

VI  =VOLT  AGE 

you  (_ptknt  )=y  1 


COMPUTATION  OF  MATRIX  ELEMENTS 

!JlT=^T07rC3-C'I17TClrC2'}' 
D12*-l-.0/(Cl  ~C2  )  /  (  C2-C3  ) 
D13=-1.0/(C2-C3)/(C3-C1) 

"l)21=  ( C2 4 C3 )/  I  C3-CU/(  C1-C2T 
022= C C3+C1 )/ ( C1-C2 )/(  C2-C3) 

02  3=  (  Cl  +  C2  17  (  C2-C3 ) /(  C3-C1) 

”  1)31  =~-C2  «CT7(  C3~C  l  )7  (  Cl-Cl  7 

D32=  ~C3  *C1/ { C 1-C2  )/ ( C2-C3  ) 

D3  3=-CI*C2/(  C2-C3 J/IC3-CI ) 


455  jD 
■A5600 
_45700 
45300 
45900 
46D00 
46100 
46200 
46300 
464  00 
46500 
46600 
'  46700 
46300 
46900 
'47000 
47100 
47200 
47300 
47400 
47500 
47600 
47700 
4  73  00 
473  2  0 


EVALUATION  OF  VOLTAGE  COEFFICIENTS  47340 

4  73  SO 

" 6  V  =1) il  ^ V  I V[j  1 2  ♦V 2  +  D  f 3 v V  3  ‘  ~  .  '  47880 

3V  -  02  l  *V  l  +  C2  2  * V  2  +  02  3* V  3  4*901 

7/  =  031*VUD32*V2  +  D33#V3  _  _  _  47920 

"  i  F  {'{'AMI  NI  ( Vl »  V2»  V  3  )  .LT  .KTOE )  •  AN"6V(  INDXETA  .EQ.4  I 5  ‘  GO  TO*  520 .  47*40 

‘FMNN.LT.Nn.OR.CNP.ir.ND  )  GO  TO  520  47530 

47960 

REFI  NED  EST  IMATE  OF  MAX  IMUM'  POWER  POINT  *  . .  47950 

48005 

.  r-  ( I NDXET  A.NE.3  5  GO  TO  395  _ _  40020 

4  I  NDXET  A=4  '  . * . . .  4804  j 

v.ALL  CRVF IT (  CRNT  »  VOLT »  PT4NT f  INDXETA 3MA  XPDW,  PMAXPOW)  48060 

t  H  AX  POW=  PMAX  P  Orf  /  CM  AX  P0a<  4  80  JO 

. .  . . . .  ‘  "48100 

VEST  FOR  ETAMAX  4  82*J.- 

4  330  0 

5  ;F (I NDXEt  A.EQ.2)  GO  TO  396  “  4846: 

IF  (  (  ET  A.  9T. ET  A1 ) «QR.{  INDXETA  .NE  .  I )  )  GO  TO  397  4 

6  1  NDXET  Ar  INDX  FT  A+l  __  __  43500 

7ALL  CRVFI’TI  CRNT,  VOtT,  PT<NT,  INOXETA70URDNET,PMAXPOW) .  *  46  7  JO 

01  FMI  N=  AMI Nl  (  ABS  (  Cl-CURONET  ),  A8S(  C2-CU:1DNET) )  *  4 8 7 4 0 

RAT  I  G -  D T  !'M  IN/  CORONET  40300 

5  F(  RATIO.LT.  0.00  3)  GO  TO  394  . 463^ 

IFdITVAR.EQ.lJ  GO  TO  410  48330 

MR11F  (6,590)  491  JO 

. TVO  fO  410 . . . ■  ----  -  ---  --  -  . . 4 920,> 

49240 

COMPUTATION  OF  CORONET  45290 

4  *  4  4  ■ 

■*  t; R'J.NP f  =  AV AX  1  {  Cl ,  C2»  0 3  )  +C INCRFG  •  ',9^.' 

•NM3*  IN-3  _  _  _  M  44 

UU  79  ■>  NY  ^  2  ,  J  NM 1 1 J  NM  3  "  45'  70 

1  -NS  ( 2  ,  NY  >  +NS  (  1 ,  -NY  )  <1  *y  •  j 

*«:*  -  NS  l  (  2 ,  NY  I  ♦  NS  U  It  NY  )  c  >  -J  3 

713  =  NS 2  ( 2 1  NY 5  <•  NS 2 (  1«  NY  5  4 


c  ■»  o  o  f.  o  o  r»  o  o  i  4>  o  o  cV 
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F 1 1  =  -i:0/(  N3  -  N  1)7  (N1-N2T 
F12=-l • 0/ I Nl -N2 ) / ( N2-N3 ) 

FI  3=-l  •  0/  (  N2  -N3  )  /  (  N3-N 1  ) 

F2 1  =  -fN2 +N3  )  *  F  I  1 
F22=-{N3+N1)*F12 
F2  3=-(N1+N2)*F13 

*  '  ~F3l=N2*N3*Fll  -  -  '  ■■  . 

F32=N3*Nl*f-i2 

F33=Ni*N2*F13 

- A=  F 11  *C  I  +  F 1 2  *C2 + F  L  3  *C3 

B=F21*C1+F22*C2+F23*C3 

C=F31*C1+F32*C2+F33*C3 

■” . .  NSRF=AMINI(N1,N2,  N3)*RF  *  ~ 

CU ROT  RY=  A*NS  RF*NS  RF+B*NSR  F+C 
399  CURDNET=AMIN1  (  CjRDNEF ,  CUROTRY  ) 


I  NTERPOLAT  ION!  FDR  NS 

10  CURDNTS=  CJRDNET  *CURCNEr 
00  42 D  NX= 1 ,  IN 
DO  423  NY=I,  JN 

. NS3  UNX  ,  NY  )  =  N$  2  (  NX  ,  Ny  K 

NS  2  (  NX  ,  NY  )  -N$  1  (  NX  ,  NY  ) 

N’S  1  (  NX ,  NY  )  =NS  (  NX  »  NY  ) 

A=D1  1*N$1(NX,  NY  )+012#NS2(  NX »  NY  )  +D 13*  N  S  3  ( N  X ,  NY  > 
0=  02  1  *NS  1  (NX,  MY  )  +  02  2*NS  2  (  NX,  NY  )  +023*  NS  3 1  N  X  ,  NY  ) 
C=D31*NSl(NX,  NY)  +  D32*NS2(  NX , N Y  ) +033* N S 3 ( N X , NY ) 
4? J  NS ( NX , NY )  =  A*C JRDnTS  +  B*CURDN  ET  +C 


INTERPOLAr ION  FOR  NAG i 

DO  433  NX  =  2 ,  I HM 1 
00  433  NY  =  2,  JUMl 

MAGI  3  l  NX,  NY  i  =  MAG  1 21  NX,  NY  j .  . 

MAGI  2 ( NX , NY )  =  M AG  I ) ( NX , MY  ) 

MAGI  l ( NX , NY  1  =  MAGI ( NX , NY  ) 

A*  01  I  *  MAG  I  I  (  NX,  NY  )  <C12*MAGI2<  N  X ,  NY  )  +0  1 3*MA  G  I  3  (  NX  ,  N  Y) 
8=  02  l  *  M  AG  1 1  (  NX,  NY  )  +C22*MAG  I  2(  NX ,  N  Y  )  +02  3*MAG  13  (  NX ,  N  Y ) 
C=  03 1 *MAGI 1 ( NX, NY  )+C32*MAGI2( NX , NY ) +033*MAG I  3 ( NX , NY) 
430  MAGHNX,  NYl  =  A*CJRCNrS+B*CURONET+C“ 

PT  KNT  =  PTKNT  +  1 
44.)  ET  AI  =  E  T  A 
GO  TO  40 


CHECK  AND  f N C R E M. f: N f  MESHXNF 


450  MESHKNT =  HESHKNf  +  I 

i 

BEGIN  RFMt-SHlNG  BY  STORING  NS  IN  DUMMY  AND  MAGI  IN  OUMAG 
IX  =  i 

00  470  Nx  =  1,  I  \  . . ~  •  '  ■ 

I Y  =  1 

DO  4  60  NY  =  I ,  J N 

DUMMY  (  I  X,  IY)  -  NS  (  N-.  ,  , 

DU  MAG  (  (X,  f  Y  )  N  !  i 

460  l  Y  = ; v ♦ 3 

4/0  I  X  = 1 X  *  3  ‘ 

COMPoFt-  OJMMY  AND  LJMAS  AT  h)+ I2D\YA.  INTERMEDIATE  POINTS 


4  v6 / 
497  4 
497» 
493  4. 
493  v: 

4  994! 

4  999: 
503V 

5  00  VI 
501  >4  ( 
5019: 
5024: 
5029: 
5034C 
5037: 
5043: 
5050! 
5060C 
5075C 
50903 
5  100C 
5)10  3 
5120  0 
513.) 

5  l  V)  J 

515) : 

516) 0 
5  ,i  7 ;  ? 
5  130* 
>  1  4  ' 
5  20'') 
5210  * 
52?  00 
52  3  00 
5  24^3 
526  0  0 
525  >:• 
527)0 
s  ?  ?  *:  :> 
W1  ; 

5  5 .  <  0 

f  j  » 

-•  t 

5)2  '/ 
33^.'. 

6  340  • 

4  $  ■} ' 

5 

5  5/.. 
5;w.v> 

6  39  '  1 

5  4  O'*  • 

6  4}  '  ' 

4  4  5 :  '• 


4  4,'  * 


i KAX=  3  *  l  NMl  +  L 
J LEFT  MX=3*JN-5 
DO  480  IX=l, IXMAX, 3 

HO  480" I  LEFT  =  I  t  I L EFf MX, 3  "  “  . 

I  R  1 GWT  =  1  LEFT  +3 

DUMMY!  IX,  I  Lfr  F  T  +1  }  =  {  DUMMY!  IX,  IL  EFT !  +DUMM  YU  X ,  I LEF  T)  +  0  UMMY  (I  X,IRIGKT 
Li!/5.3 

:/JMAG(  IX,  ILEFT  +  1  )  =  !  DJMAG!  IX,  IL  EFT  )  +  DUM  AG!  I  X  ,  I L  EFT )  +DUKAS  II  X,IRIGHT 
1 i /3.0 

DO  MMY  (  IX,  I L  E  F  *"  +  2  i  =  (  DUMMY!  IX,  IL  EF  T)  +DUMM  Y(  I  X  ,1 R  IGH  TU  D  UMMY II  X ,  I  RI GH 

'ni/3.3 

-  DO  MAG  (  I  X  ,  I  l.FFT+2  )  = !  DO  MAGI  IX,  ILEFT  I  +  DUMAG!  l  X ,  IR  IGH  T)  +  DUMAG  ( I  X, !  RI  GH 
IT  })  1 3 . 3 

COMPOTE  DUMMY  AN  DOM  AG  AT  VERTICAL  INTERMEDIATE  POINTS 

i  l.OWMAX=  3  *  IN -5 
l  YMAX  =  3  X-JNMl  +  i 
DO  MO  I  LOW- I , I LOW MAX • 3 
DO  490  IY= l, IYMAX 
IHiGM* I  LOW +  3 

DUMMY!  1LOW+1,  !Y  1  =  {  OjHMYC  ILDW*  IY  I  +DUMMY  (  !LOW,  I  Y )+D UMMY (I Hi  GH  ,1  Y)  ) /3 
loO 

DUMMY!  ILOK+2,  IY)  =  !DUMMY(  1LOW,  IY  }  +DUMMY l  IHI  GH,  I  Y )  +OUMMY  ( I  HI  CH. ;  1  Y  ))/ 
13.0 

DU  MAG l  I  LOW  +1 ,  IY  )  =  (  DJMAG!  ILDW,  IY  )  +DIJMAG!  ILDW,  1  Y  )+DUMAG  (IHIGH  ,1  YU/3 

1.0 

J  DU  MAG  f  l  LOW  +2,  IY  1  =  f  DUMAG!  I  LOW,  IY  I+OUMAG!  IHIGH,  I  Y  )  +DUMAG 1 1  HI  GH  ,  i  VI  )  / 

iO,  0 

IRANbFER  DUMMY  l M f 0  NEW  NS  MESH 

)  N-  3  ❖!  N-4 

JN=  3  *  J  N-4 . 

’■*0  5  30  NX=  1,  IN 
DO  500  NY  = 1 »  J  N 
IS  !  NX  ,  NY  )  =  OJ  MMY  !  NX  + 1,  NY  + 1  ) 

I RAMSFF.  P  DJMAG  INTO  NEW  MAGI  MESH 


f  H-3  *  I H-4 
UC  510  NX  =  I ,  in 
00  510  NY  =  l ,  J  ••! 

■•'Aol  <  NX  ,  NY  3  =  DJMAG!  NX  +2,  NV  i 

RESET  l  ANC  J  10  f  Hi  NJMEER  OF  INCT  C  M£  N  T  S  IN  THE  NEW  MESH 


i  -  1  +■  I  +  I 
J-  J+  J+  J 

UT.GIN  NFW  MF.N.H 
GO  TO  30 


SUMMARY  WRIT  f  UJf  OF  VOLTAGE,  LL-RRl  \  T  AND  EFFICIENCY 


•' .  AM  -  Pi  MO  -i 
MAT!  (A,  500  ! 
i)r  ■> :  0  1 MDX  -  J  »  PT  \  IT 
YOU  T  mp=v?;lt  (  INOX  ) 
IDT  MP~C.k\T  (  1  No*:  | 


55100 
55200 
553D0 
55400 
55500 
55500 
55700 
55800 
55900 
5  5  J  0  0 
56100 
55200 
56003 
55400 
55500 
56500 
56700 
56800 
56900 
57000 
57100 
5  7200 
57300 
57400 
57500 

5  760  0 
57700 

6  73  0  0 

)  %  *  \f 
>  \  0 "/  .* 
58100 
5  8200 
583)3 

5  64.)"/ 
5850) 
56600 
50703 
5B800 
509)0 
59j0) 
69100 
69203 
59300 

k«-it  l) 
3  9  •> .,  J 
59703 
59.303 
5993  j 
6000 
6010  ) 
602  3  v 
60)0  , 

6  04  0  * 

*■  j  0  '  1 

6  :■  : 


w  ') 

637  >  ’ 

4  '•  V. 

61 

513). 


_-8n.-  __ 

tr  A=  CJ  RUT  MP**/  OLTTMP/PQwMAX 

ifucjrve.eq.l)  go  ro  530 

WRITE  (6,610)  INDX, C'JRDTMP, VDLTTMP,  ETA 

'  •>  1  GO  NT  I  Nj  E 


"UMMARY  OF  OJTPJT  RESULTS 


JCC  =■  \>  CLT  ( I  ) 

VFsVOC/VGAP 

CL'ROSC=  (-BV-SCRH  3V*BV-4.O’0‘AV*CV )  )/2lO'/A'V  . . " 

!'  T  AC  01.  =  C  J  RDS  C/  C  J  R  DM  AW 
-F  =  P  MAX  POW/  V  0 C/  C  J  »  DS  C 

rr  ACELL=  PM  AX  PCM  /  POWMAX  ’  ~  • 

WRIT  U  6,6  20)  v'OCi  VF,  CURDSC,  ETACOL  ,  5M  AXPOW,  CMA  XPO  W  ,  VMA  XPOW,CF  , 
.FTACELL»XO,LN,LP,LI,£NR 

FORMAT  (/,  IX,  SHy'OC  =  ,  F6.4, /,  IX,  5HVF  V  ,  F6. 4 , /,  1 X,  9H;  URD  SC  =  T 
u:;o.3,/ ,ix,9Htr acol  =  ,  F6.4, /,  ix,  iohjmaxpow  =  ,eio.3»/, 

1  X  ,  1  OHCMAX  POW  =  ,  F10.3, /,  IX,  IOHVMAXPOW  =  , F6.4 , / , 1 X, 5MCF  =  , 
if  6.4,/  ,  1 X  ,  10  HET  AC  FL1  =  ,  F6.4, /,  IX,  5HX3  =  ,E  10.  3,  /,  IX  ,  5HLN  =  ,610., 3 
't  ■>/  ,1X,5HLP  =  ,  EI0.3,/,  IX,  5HLI  =  ,  E10.  3, /,  IX,  6HENR  =  ,E10.3) 

SI  ORE  OUTPUTS  IY  VECTORS 

-  OCD(  I  0)  =  VCC 
VFl){  I  0 )  =V  F 
:.UftOSCD(  I D ) -  C J R OS  C 
w’J  MO  MX  D  (  !  0 )  =  C  J  ROM  AX 
4  T  AC  OLD  I  101=  tr  ACOL 
"‘YPOWOt  I  U !  =  PMA<  POW 
C:-'>OQW,.)(  1 1)1  =  oiy 
vf'XPQrflH  l  0 )  =  V  V  A<  P  OW 
Cf  LI  [  0  j  -  C F- 
cl  AO',  i  Di-l-T  ACPLL 
OLNt  ID)" ML') 

JN!.f*U  0)  =  \LP 

,  ;<ouo)  =  m/o 
IF i I  3  PTC .EC. L  )  YJ  T 

S r  IECT  Ntrl  M fc  S  ri  p  AR  AM  ET  £R h 


730 


C  . 


If  (  iC.cOO  ) 
JFUC.EC.2) 
i  i  u «;.  t.  s. 5 1 
lr-  (!  C.  EC.  A  I 

:  f  1 1  c.  t  ;,S! 
irnc.Qc.s) 


VL M  L  =  MLMT~  I 
L  M  )  -  N  L  N  T  +  l 
)«.  P 1  -  M  pr-  i 
v,i.  0  !  -  P  T  ■<  ’ 
’KOUMXf.T-  I 
NXf.l'NXOn  I 


)>'<VE  OJTFJTS  CEr":  tV  At  J  A  T  Ft)  f  OR  NEW  ME SH  PARAME Tb R S 

IEI*  IEVAL  (NL-'l,  Ml  PI,  MXJl> 
fF(in.EG.P)  GO  10  ‘> 

'S  POWER  OFJTPJi  F'lR  MfvTST  UrGt«  GREATER  THAN  MAXIMUM  "OKER 
U’TtVl  TRY-  prtWIt.'JT,  MbSMF., 


61400 
61460 
61543 
51S03 
61700 
61900 
61.3  30 
62030 
62100 
62200 
62300 
62400 
6  2500 
62600 
62703 
6  23  0.' 
62300 
6  30  00 
63100 
63200 
63300 
6  34  00 
635  00 
6  ■>  5  0 
6  3700 
6  3300 
6  33  MO 
6  4000 
6  4 100 
6  4  W  U  J 
643  00 
64400 
6  4  6  0  0 
h  4  *•  0 
64700 
64  530 
6  4-00 
6  50  GO 

6  ‘3 1  .or 

O  ^  u  'V 
!> ') V' 
66*  !r 
65..  O 
t*  ■> 6  0 . 
6  57  00 
649'. 
663::: 
>?oo. 

V)>  3. 
647  ,K 

5 

6  .>  4  j  . 

4  45.1 


f  _  * 


I 


c  C  V’PIJ  GO  TJ  720 

'  .  GO;  V-:'  i  G  J  :  Pv  \ 


(  J  ■» 


1  U«>  1  O  It  V 


<  i 


«s 


*.  •; 1 
d  o  • 

^  i  >  • 


-85- 


NXOT  =  MX01 
IDXXOM=3 

_ I T  =  I C 

I3PTS  =  0 
GC  ro  700 

;  '  HAS  THE  POWER  OJT  PUT  BEEN  .EVALUATED  FDR  THE  -SIX  ME  SR 

;  ARRANGEMENTS  ADJ  ASCENT  TD  THE  MESH  MAXIMUM  POWER  POINT 

* 

?2  0  fCMCVl  . 

IFIIC.EC.7)  IC=1 
NLNi  =  ML  NT 
NLPL  =  NLPT 
NX  OL=  NX  OT 

IFUC.NE.in_  GO  T0_  700 
C  EVALUATION  OF  OJTPUT  PARAMETERS  FOR  3  ADDITIONAL  POINTS 
13  PTS  =  I 

JLNP=  IEVAKNLNT  +  l.NLPT,  NXOT  ) 

JLNM= I  EV AL ( NL  NT- 1 , Nl PT, NXOT  ) 

~  . . JLPP=IEVAL<NI.NT,NLPT+-1,NX0T  1  . 

JLPM= I EV  AL { Nl NT, NIPT- 1, NXOT  ) 

JX  0P= I EV  AL I NL  NT , NL  PT, NXOT  + 1  ) 

J  X  OM  = I EVAL( NLNT , NLPT, NXDT-l) 

I C  LN= -I 

IF  (PMXPCWDI  JLNP)  .GT  .P.MXPOWDI  JLNM  )  )  ICLN  =  I 

. RICLN=  ICLN 

I C  L  P=  —  1 

IF  (PMXPCWDI  JLPPI  .GT  .PMXPOWDI  JLPM  )  )  ICLP=1 

Ri  CLP=  !  CLP 

ICXOs-1 

IF (PMXPOWO< JXOP) .GT.P^XPOWDI JXOM  )  )  ICXO=I 
RI  :xo=  ICXO 

730  LNLP=  I  EV  AL  ( NLNT  +  I  CI.N,  NLPT  +  ICLP,  NXOT  ) 

INXO-  I  EV  AL  (  Nl  NT  +■  I  CL N,  NL  P  T ,  NXDT  +  ICXD  ) 

“  LPXO=IEVAL(NLNT, NLPT  +  ICL P, NXDT  +  ICXD  ) 

Nl. Nl  =  NLNT  +  I  CL  N 
NL  PI  =  NL  PT  + 1  CL  P 
NX  01  =  NX  GT 

I  F  H  NLP.  EC. 0  )  GO  TO  e.i 

1 F  (  P  MX  PfW  D  (  T  F  V  )  .  L T  .  PM X POW  DI  L  ML P  )  1  GD  TO  710 

NL  PI  =  NL  PT 

NX  0!  -  w/f'T  T  C  X  n 

1KNX01  .EO.NXOMIN)  GO  TD  7 32 
IF  (LNXO.FO.D  !  GO  TO  5 

IF  ( PMX  PfIW  D{  I  FT  } .  LT  .  PMXPOW  Dl  L  MXO  >  1  GD  TO  710 
Nl  Nl  =  Nl.  NT  > 

NL  PI  =  NLPT*  I  CLP 

IF  {  LPXO.  EU.D  )  GO  TO  5 

IF  (PMXPCrfDI  J  t'  ).LI  .PMXPOWDILPXO  ))  GO  T-J  710 

C  ■ 

f  EVALUATION  OF  COEFFICIENTS  FDR  POWER  IN  SECOND  ORDER  EXPANSION 

CALI  CCUf 


t».  w'OD 
6  1 5  3 
6/533 
6  7  3  j  J 
67633 
6  73  33 
6? '3D 
67333 
6^03 

5  A  3  3  3 
681  33 
08  233 
68  SCO 

6  8  '■  33 
6  C  5  3  3 
6  H  ■'  3  J 
6  3  7  J  0 
6  8333 
68933 
66033 
6  ->  i  33 

5  9  2  3  3 

6  93  30 

5  o  no 

6  3503 
66833 

5  9  >  j  0 

6  9.>.  k1 
6  U33 


i  • )  <,  1 1 J 

7  (•  ■■  \) 
7  0  i 3  7 
705.0 

706  33 

707  33 
7  0*<,>3 
7 1.  -1  ^  3 
7  ; ...  3  3 
7  i  l  >J 
<  i  r-  3  3 
(1**3 
7 ;  i  o 
r  i  -  3 
7  5  :•  3  3 
7  1533 
?  1  <\‘3 
7  1  3 

71  i  3  3 

72  3  * ) 

*  y  l  v* 

7,  ?33 
7/  .'P 


JSc  OF  l  I  8  A  A  K  v  s  J  t'  -l  i  VI T  J  \  t  T  1  (  Oil)  l  N  h  TPJl  A  TIO  MAXIM')'-. 

POOTR  POINT  ’  7. 


S» 

/63  CALI  Sil-roH  AM,P,M,Xw,  3,  UVX  h  ,  »,  f s;.  , 
*Fl  I  6  l  .NO.  01  GO  TO  7r.Q 


/ , 


">  O  O 


t  .. 


65 


737 


4  0 


' » . 


..  o 


»  r  v  I  DX  X  OK.  EQ  .  2  )'  X Ml  3 ,  U  =  6 

0  7  35  NX= 1 1  3 

XMlNX,  1  )  =  -0.5*XM(  NX,  1  5 


EVALUATION  OF  OUTPUTS  AT  INTERPOLATED  MAXIMUM  POWER  POINT 

r'*MAX  POW-0(  PMX  POrf’Oi'I  FT  ),  BM*  AM,  XM  ) 

NL  N=  NL  NT 
NL  P=  NL  PT 

NX  0=  NX  OT  *  '  '  - .  . 

NLN=  NLN< X M ( 1  ,  ]  } 

NLP=  NLP+XMI? ,  1  } 

NXO=  NX0+XMI3 , 1 )  . ~~ . .  ~  - .  ‘  ’ 

:rUTfc(6,6&05  NLN,NLP,  NXO,  AM  ,  BM 
FORMAT ( 3FI3. 3 ) 

!  r-  (  a mci  ,  2 ) . e o . ovr  go  * r o"7 3 r  - - 

V.ALL  LIGENI  AM,  El  OVAL,  3,  IDXXOM,  J. } 

.'RITE(6,560) 

NR  I T 6  (  6 , 660  5  E IGV  AL, AM  “  ----- 

CALL  COEFFICMXPCUD) 

CM AX  PON  =  0(  CMX  POW  Dl  l ET ) , BM, AM, XM  ) 

-mLCCDEFFIVOCD)  . .  *■  . 

VOC=  D(  VOCOI  I  ED,  BM,  AM,  XM  ) 

call  ccef-fi curds cdi 

omosc=  D(  CURDS  CD(  IfcT  ),  BM,  AM ,  X  M  I 
'  ML  C(JEFF(  CU  RDMX  0  5 
N'0MAX-0(CJRLMXD(  I  FT  ) ,  BM  ,  AM,  XM  5 
v  -=  /OC/VGAP  '  '  "  ‘  ~ ' 

,  .  /‘.OL=CJRCSC/CURCMAX 
’  .'AXPnu=PMAv  pin;  /  CM  AX  POd 
’'MAXPU.V/v  OC/CURDSC 
'MU  -VfoCF'ET  ACOL 
*  <  =  VGAP-L  I  RDM  AX 
'  ,.‘i#K<2 .0 

m  pm;*,  .u 

D-NXO*!! 


.<  -> 


Hours 


Ml  »'(«>, f.2«>J  U,  V,  L  I,  CUROOPT,  ISP 

OR!  CM  .  /  ,  /  ,  i  x  ,  4  hH  ,E9.3 ,2X,4tM  -  *F  V*  3,  2X,  5HLI  =  ,E9.3,2X, 
ou;uRIJLr>I  -  ,  f  •>  .  3,  2X  *  6H1SE  =  ,  II,/,/} 

M  ;  6  ,  A  it) ) 

■iR  *Ar  (  1  <t  .Dili)  ,  UHLN  ,3HLP  ,  3hXJ  ,  6HVUC  ,  fcOVH  ,  1  OHC  USD  SC. 
.'.AC  ,  *()Ur  ;  ROM  AX  ,  JOMCMAXPOW  ,  CHVMXP  W  »L0IPP4XP'JW  , 

OCT  .  5  .{PT  M.F,  /  1 
I  :  7-3  IP- I  ,  ID 

■.O'  !  I  i  {  s  I  >  '•  U !  IP,  JNL  it  JP  >,  JNL  -M  IP  )  i  j  M  a  T  (  IP  }  ,  V  >C  0  (  I  P  I  ,  VH  0  1 1  P  J  , 

:  i  T'SCDl  I  Pi ,  I  )  ACOL  Ijl  f  P  ) ,  C  IR  OMXDl  IP  } ,  CMXPOWD  IIP),  VMXPU  '*0  1 1  P I  , 
O-'vPi/;,'.  1  P>  ,  ,  S-  J{  ,  L  f  AO(  t'A  1 


i  A 1 3 . "  i  am.,  rn; 
I  i  i*  (  ,  *>  5  0  J  I  .  ,  L  I ' ; 

rX  t'Urt'  i  P M  * x  i>'  ■,{ , 

M'Af  (  /  ,  /  ,  r< , 1  IN  M 
7 :  v 
;  ,  /  I  "  '  • .  '■  - 

/  ,  !  '  •  *  *  i ' 


»3f  f  6.4,  2EIC1.3,  rf./,,£l  ;-.3 ,  ?r-6.4) 

40,  7  JC,  v I  *  Cu><DSC,C  :  K  .CfROMA \,C  VA  XPOW, 

AV  ,  '  c, rr A  FLL 

;u.  3,  / ,  6i*  !.  r  l:  I C  •  3,  /  ,  5H  xo  !  L3  ,/  , 

r  ;!:N,  C  -  .3, /,<>-.  P  ULU  =, 

•  •  10  .  T.  /  (  I  ‘  .-  =  . .  ■  .  X  ,  i  ,  j  i ;  .  V  .'•*»<  PC  .•< 

•  no :■  \  .-to,  -,t  6.4,/ 


‘  i  r .  ■ 


f 


» 


7291  , 
730); 

7  31  or 
7323; 
7  333. 
7340; 
7350; 
7353: 
7  3700 
7  38  OC 
73900 
7403  0 
74100 
74200 
7  43  00 
7  435  0 
74400 
74500 
7  4  5  j  0 
747  03 
74330 
74500 
7  30  Go 
7513  J 
7  32  3'. 
753.0 
7  54  00 
7  5  3  .)  J 
756)0 
7  5  7  0  0 
7  53  M 
7  5  n  :■ 

750  0  4 
761  Jo 
7b?  j"» 
7  5)5 : 
7  64  0.1 
7  65  00 
7  66,-' 
767  33 
7  63  0  0 
7  63  >0 

7  70'' 

7  7):  ' 
V  '( (  >  •• 

i  l  'j  V/  .< 

77  40  • 

v  -  •»  ^ 

t  !  -J  I  > 

I  i  5  -*  J 
7  7?)) 
?  I  .00 
7  7  )  0  ) 
7000  5 
7'P  3 


■  ' , 

..T  (/  ,  ! 


'■  ,•  t 


I  -  f 


i  - I 
7  •  . 
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AM.',  11  =  0  786 30 

^(3fl)  =  AM(I,  3  )  =  0  735%j) 

*  >v  1 2  ,  3  )  =  AM  (  )  ,  2  )  =  0  73580 

«’ *'  TO  733  ‘  '  "  '  78550 

787  00 

(!X,1  7HS  iCCESSFJLL  DX  IT  ,/}  78300 

.'M.'AT  (  IX, 5  (  2X,  H  3.6  )  )  7  9900 

'  0K.vAT  !  1 X  1  7  0000 

!  f'WMAI  l  l  X  ,9  Hwt  S  IKNT  = ,  i  4,  9H.  OO^KM  I  =,IA,/}  79100 

'>■  ■JM’-'Ar  !  /  ,  /  ,  I  X  ,  9  1  =  ,  12,  IX  ,  OHP  A  ThK'-i  7  =  ,I2./>  *  *  '  79200 

l  /  ,  /  ,  1  X  ,  ?  7  H  ^AXIWJM  r  FF  ! C  I  LMC  Y  VALUES  ,/,/)  7930.) 

0(mr  (/,/,  l*,5HMK\lT,  6X ,  ?HUH0MET,8Xf  7HVQL  TAGE  ,  8X,  3HE  T  A  *  / )  7)400 

F M R  M  A  T  (IX, 15,4 ( 9X,  1  10.31  )  795  00 

f  C8MAF  (  i?X  ,  f:  1  2.4  )  79550 

~ 1  •  u  7  c>  5  0  0 


-bb- 

jui.huJi  INC  LOCK- (DM) 

COMMON  I  OT  ,  JL  MPt  J l. NM v  J L 0 P f  JL  PM ,  JXDP , JXPM ,  I  *.'LP  ,LNXO,  1  PXO ,KI CLD , 
-  1  L l  P t  P I  C <  C .  BM(  3 »  I  )»  AMI  3,  3  )r  J  'U.M  (  50  ) ,  JN|_P  I  50  ,  J'JXOI  50)  .ID 

DI  MENS  if:.  DM {5J)  .  -  - 

AM  I  i  ,  1  )=  (  DM (  JLNP  )  +D*M  JLMW  )  }/2.0-()M{  [F  r  ) 

>iv  (1  ilIM  DM{  Jl  NP  )  -0M(  J  L  N  M  )  1/2.0 

*lM<?  ,2  )=  (  ()M{  JL^O)  4-DM{  JiPM  )  )/2.0-tJM(  I  £  T  ) 

•  M  ( ?  »  1  )  =  (  UM(  J  L  PP  )  -C  M(  J  •_  P  M  )  )  /  2 . 0 

A*  (3 ,3  )  =  (  0M(  JX  )P  )  +  Qv  (  J  >:■ )  y  )  )  /  2 . 0-  Ov  {  It  T  ) 

r.M  M  , 1  1=  (0M{  JU3P  )-0V(  J  (DM  )  )/2.0 

A  M  I  #  ?  )  -  AM  (2  *  l )  =  (  CM  {  L  \’L  P  )  -  DM  I  1ST  )-BM  {  1 ,  1  ICLN-6/K  ?  ,  1 )  r  IP 

)  -AMI  I »  i  )  -AM{  2  i  2  )  )  /  2 .0/R  I  Cl\>/^  I  CLP 
AM  (  L  »  M  =  AM(3,  1  )-3. 

AM  (2  ,  l  )  =  AM!3f 2)  =  0. 

IT  (  (LMC.K.D  )  .0K.{  LPXD.f.3.0)  )  C  T  J  ■{  v 

\"(I  ,  3  )  =  A  {  3,1)=!  iM(C\<0  1 1 T  )-3M(  [,  )  ICCM-:5M(3  ,  L)*fUCXC 

•  M)  “AMI  3,  3  !  ) / 2  .0/9  I  CL  M/9  IC  XM 
"r  (?  *3  )  -  AM (  3  j  2  )  =  (  Dv  ( L  P< 0  )  -DM(  I  FT  )-SM(  2*  I ) *9  ! CLP-3 M <  3  , 1)  * P.J  C  XO 

'  AW(2.?I-AP(  3,3)  )/2  .O/RICi.P/3  iCXJ 
‘P  T  JRN 
‘  U> 


797,90 

7/933 

7  9  9  J  ) 
bOj  )  _) 
Bui  30 
P0239 
BO  3  2,3 
f'CK  3  3 
‘’05  2:3 
9053  ) 
AO? ' 
50?  3  3 
80520 

8  05  A  3 
BO?  5J 
fj  O  :  3 
8)039 
fi  l  i  33 
8)293 
8’  j  02 
81400 


-89- 


r 


FUNCTION  C( diet, »m, am.xm j 
01  MENS  I  CN  AMI  3,  J  ),  TiM{  3,  LI,  XM{  3,  1 )  '  ‘U  '  IQ 

0=01  (IT  'U  JO 

00  8  30  1-1,3  -  -  ......  ‘3  mo 

OsO»BM(ItH*XM(  !,!l  n]!'  'J 

00  300  J=  1 ,  3  «1/J0 

D*D+ AM(  I  ,  J  M(  I  ,  1  )  *XM(  J,  1  ) 

RETURN  >i/J  3*/ 

end  0 

R2.I00 


-90- 

PUNCTHJN  iEv'ALl  -VLNl.NLPl,  NXOl  )  H, 

CCRMCJN  IET, JINP,  J  L  P  P,  J  LPM  ,  JX3P  ,  J  XQM  ,  LN‘  ?  ,LNXO , LP XO , R I  C LN ,  8 . 

lRICtP,RICXC,J3M(  3,  1),AM<  3,  3),J'JLN(j>0),J\'LP(  50,  ,  JNX0I50)  ,ID  RTS* 

I  p=o  -  -  .  .  '  '  a.  / 

i  e  v  a  i. = o  R?a;: 

810  IFUP.EC.IC)  KErjRN  8290 

I  P=  I P+ 1  '  "  ’  8  VjJ 

IF  (NLM1  .NE.JNLNf  IP)  )  GO  TO  810  831'' 

IF  INLPI  .NE.JMLPI  IP)  )  CO  TO  8  10  _  8  32J 

IFINXOI.NE.JNKOI IP))  GO  TO  810  '  “  *'  8330 

IEVAL=IP  8  J  A  . 

RETURN  8  3'.;  j 

END  8  360 


-9i- 

f  i  'it  CK^IK  CRNl  ♦  V  JL  r,  0  IKNT,  INOXE  T4,CIH0N6  T  tPMA  XPOW) 

’  '  rt  ,F.  R  PT  KNT 

■i  t.  N  >  I  CN  VQL  1(33  ),  CRNT(  30  ),  C'U  30,  30)  »  Tt  3C)  ,  FIM1  C30)  , POWER  HC) 

-  \  .  (  JO  } 

N  >  PCwtR  I  N  T  ERMS  G F  CURRENT  AN 0  VOLTAGE 

■'  1  i  “  i  ,  P  T  KNI 

K(  l  )  =  V  C!  r  (  1  )  -CRN!  (  I  ) 

ltd  PCLVNCWUL  EXPRESS  I  )N  I-3R  POWER  in  terms  of  skewed  CURRENTS 

K|  ■  --  AKAX1  (ONi  (  n  T  K  NT  ) )  L  5  N  T  (  3  T  ^  N  T  -  1  ) ,  C.  R  N  T  (  P  TK  N  T-  2  )  ) 

■  .  '  A M A y  1  (  P')„  F'U  p T < N !  ),  PjWf  R(  PTKN  r-  1)  ) 

/  >  I-!,PT^NT 

■'.1,11-1 

•*.,  (  1  )  :  C  R\T  (  ll-SNf  «(  t  ) 

i’-'  j=;.’,prKNi 
•'  (  1  ,J|-CMII,  J-l  )*Ct<NSI  !  I 
C  A  v  L  L  I  M  Oil  CM,  P  )WT  Pi  T»  JO,  PUNT,  1,  1LE  I 
■  7  !  i  r  ( 'i ,  3  J  III. 

<  <*AT f  S  H I L  T  =  ,  I  I) 

f  I  NO  THE  MAXIPJM  POWER  POINT 

*)'  .  I=/>,PTKNT 

*■  "I  !  -1 

i  vi  <  i » - n  n  *w  i w ) 

•  or  i  --  ■.  -ns  (pi  ^ n r  i 

■  M  ,  -  C  <’Y,  !  "[  KM  -  l  ) 

•»  I  -  i ,  3 

>■  1-  -•  ‘.MM-lPTNOr  i  i  I-U,  CRNT  1  ) 

'  -  -  I  >-  i  !  M  MIT  ,  !  l»*  I ,  CRN  I  ) 

•  .  (P-'K  (  .  I"’  )  >> !  I ‘i  A 

.  .I  F  '  PPJ  <  (.*■  Nl  ?-Pp  lJ/(  I’P  7  J 

!■  \I  (  R\T  l 
'  !'M  I  ■  t.UPO'jH 
.  NT  IN)' 


-  ;  '•}»  JT  At  l  I.N  'J r  MAM1  »*'  ?  ■  »M 

P«  r.  -  A  , 

.‘Ml  I 

•  K  1  J-5  ,r-T>.NF 

WPAXP’lrti  'R’A.  .Mo  till  )  c ,,  j  ; 

.  »>!  j  •* 1 *«  j  •  ,1  r 

iAK!  W  t.JWFFNT  AT  .,jnl  i  r‘  ) !  .  F 

T  •  C  1 ' •’’’N  i  T  *5  \  f  n  «!*'’*  A  •  ‘ 

;  i  kn 


8  3  /  J  j 
El  3  n,o 

R  i )  ■> ; 
84000 
3410;) 
84^00 

84 3  3  ; 

844  DC 
84  5  00 
841,00 
34M0 

84  3  00 
8  4  000 
85000 
8  o  1  O’, 
B5/30 
8  4  3  00 

85  4  TO 
355  00 
854  30 
857  0  0 
853  G, 

83  no 

4  f  )  ; 
54  i  J  » 
8  i\  ?  .*  1 

Pm11 
8.4  .  ,  ) 
fj(i  >  5 
8  45  j  / 
8  4  7  0  0 
VjH.i  ' 
86900 
*  r  .  1 

V 1  0  J 
■Mr  ■  1 

f, :  - 1 > 

U  f  ^  *  ’  ‘ 

‘  7  5  : 
3  i  r,  .  > 
k  /  ;  ;  • 
P  /ft  "i  ' 
w79  ) ; 

/•  i'.i  * 

F  ,i  0  t 
3  ■*.  ?  0  .) 
"  ►«  X  0 .) 
>'  J '  0  1 

"55  •> 

3  .ft  ,  i 

8.-.r  \> 

H  35  j  J 


„  .>..*«  '(■  ''  1  .111  r'\ NNT  •  1  l  v.  1,  LJ-iONfr  t  ) 

v'Fmtr  prKN'r 
'V.  Siir.  T  IKK  30  } 

’V  =  » 

;  -  I  f/EsO, 

Is^l’TKM 

y,._  ?>,»,?  JM,  tT  !  ”  1  (  !  }  *-C  Iy  ? 
t  ,J?  r  C  t  'V  #CJ  !'  ONt  r 

r  )  9 


fib  )  >1 

0900 
0  >  ]  » 
0  9.'  j 
89  5  3 
99'.;) 
0  959 
b  9  ; '.) 
0970 
090  5 


l 


l 


(  /  .OOF-O  j 
' ,  .  r>r)i-D) 

'  /SJDF-J} 
>.  .■ OOOE-Oi 
s  ,-s'ioe-)3 
"  •  0  L‘  -  J  3 
.  7 >00 £-03 
l  t.')0F-02 
1.  I  jt.OE-32 

1 

u  ' 's?e~a2 

.  . siOOC-02 
; .  r>  .>  n  t- -r»  a 
t .  7  6  JO  L—J? 

1. nl 50G-02 
2. 0090C—02 

!  Jr>01 — 0  2 

2.  w  30F-02 
?.  ‘>ooo  e~02 

'  .  '>240  t-3? 

2  *  />  90  1-0  2 
.  >?  Of-O? 

*. .  000!  -02 
'  2'. Of-O? 

^  DO  fir-3? 

\  ’50f  02 
?>/jor-02 
1  -  ’.2  5  9  f -02 
•  >.)0 1-02 
/OOf-02 

■> .  >  >oo  r--0  2 

.  /S')  f  -0  ? 

.  ’  ■  O  f — 0  ? 

.  .  W '.Of— 0  2 
s  DOE-02 

•  >  -JC-0? 

'  . '  >  m-02 

<>  l .  5)1-0  2 
•-  ‘  » .  ’  0  f  -  J  2 
...  I  /\>0f  -J? 

*» .  2>./.f  f-J2 
5.  3  7501  -02 
S,  SQOOf:-.)? 
s  st  so 

•  .  ■  >00  f-  -02 

*7 ‘Of— 02 
S. SU00F-J2 
s. >  2 so  H-0? 
_  0)1-02 

•>  ■  Of  -.i «’ 

,  >?<.  •  » 

*  *  '  t  I  .  ’  ✓ 

s,  *  .»  0  f;  -  >■  2 

;-n-)2 

s.“ /SO  f— 02 

0  t 


1 • 0  6S  4  ?:  1 
0.31  i  3  G  1  9 
3.4637G19 
3.777SH+19 
4.0326F+  L  9 
4 . 2 1  5  4  f  4- 1  9 
4.  35 )?£-  19 
4.'*7?:r;+L9 
4 .  S  7  6  4  f  +  1  9 
4  .  64  2  4  E  +  1 9 
4  .  7  3  f )  4  L  +  1  ) 

4 .  no  6  3  f  +  1  9 
4. 86  7  fl  Ft )  9 
4.02 3  Oft  19 

4  .  )  7  S  6  [  ♦  I  9 
5 . 0  2  3  4  4  4-1  <> 
5.96  8  1(_  4  I  9 

5.  u  )  7  (-4  19 

5  .  1 4  6  8  F.  +  1  9 
5 . 1 9  5  6  F  1 1 9 
5.2202F  +  1  •> 
5. 253  1  ft  19 
5.28  4  2  (.  +  !<) 
5 . 3 1  3  7  f  t  J  9 
5.  34  18ft  1  7 
5 . 3686f;  +  l  v 
5 . 35  4  2  f  +  W 
5 . 4  i  3  7  f  1 1  9 
5.44211119 
S  .  46  4  6f.  +  1 9 
5.4862C+1  0 
5 .  SO  S  J  f  ♦  1 V 
S. 5? 6 9  ft  19 
5. 54  SI  ft  19 
S .  54  4  6  F  t  j  9 
5  .  S  3  2  S  F  1 1  / 
5.5V  t  H  f  r  1 9 
5.41 sSf +  19 
5. 6  3  2  ?  r  j  •; 
5.64’)  j.  t  J  9 
■>  .  6  S  t  S  f  t  |  9 
5.67» 211  +  1.9 

2  4  H  <  ;  9 

5,79'.  ?C<  19 

S  .  7 1  ft  1 1  j  i 
S  .  7  3  •>  8  r  *  1  9 
S  .  7  4  “>  •>  L.  +  1 
s .  r > ;  •)  1 1 }  v 

5 . 77  )  )  t  .  1 9 
5.74  1  7ft  |c, 

5. 70 3  2ft  1  9 
>»S  >v  )ft  19 
5 . 4 !  ■>  ?  M  l  9 

r  .  2  2  ‘  H 

-  .)  <  ,  _»  -  •  . 

.  tit 

S  l|  +  * 

>  ,  6  4  S  j 1  t 
s  .  *5  f  4  -  .  ;  > 


2 . 6Q68G  +2 1 
1  .674  5F  +  21 
l.?834ft2 l 
1  .0335fi  +  2l 
8 .60 96G  « 20 
7 . 3683ft? 0 
4 , 4  4 1 9  f +20 
5 . 7?68f +20 
5. 15  791-4  20 
4.49355+20 
4.3059ft  2  ) 
3.9744ft 2  o 
3.492  It  +  20 
3 fc4 4 3 «f +20 
3.224  7ft 20 
3.02961+20 
2. 8 54 7f +20 
2. 6970ft 20 
2 . 554  Of  +  20 
2.42  371 +  2  0 
2 . 304 6f  <20 
2. 1952G+20 
2.0945F+20 
2.0016+  +?-) 

1  , 9 1 SSf +20 
1  .8  3561 +20 
1  .  76  T'f  +20 

1 .691 9F+20 

1 ,6?7?r  +20 
i  .5666c  4  20 
l . 509  76*20 
1  .4  5641 +?o 
l  .406  1  f  +2.0 
1 .358  88 +  20 
1  .  3  i  4 2 f +20 
l  .2  720*- +  20 
l  .232)  t  *  0 
1  . 194  38  +  20 
l  «  l  5  H S  t */u 
1  .  1  2  4  4  f  •  2' 

1  .'■'92  1+  *  20 
1  .061  \‘  -2  ) 

1 .0  3?0i  +20 

1  . 0 0 1  )t  •  1 0 

v  »  V  7  i  0 1  t  1  : 

9 .  s  :,j2t  +to 
74* 4  +  i  9 
1.041  31  +  i  4 
8.83  n  Of  +19 
6.604  f  ♦  l  4 

3  »  I  *  ^ 

n  .  < !  2  •  i 

I.ft'.j/!  •  !  ; 

/  -./  >->  if  , )  J 

/  .48  7 >•*  '  1  9 


co  cv 


7.  25^0 L-02 
>  , ^7^0  E-02 
■■  .60006-32 
7.32536-32 
? »  7530E-32 

7 .  ^1750  f£-0  2 
8.00006-02 

- 1 250  e-D  2 
.2530  6-3? 
E. .  37 50  f-02 

H .  5033£-32 
B.5?f>3t-32 

8.  7-330E-J2 
li .  3  750  6-32 
9. 00006-32 

9.  L 250  £-0 2 
9 . 2^30  C-J2 
9.  3  7*306-32 
9. 50306-02 
9  .  fj 250  6-3  2 
9. 7500 6-32 
9.37536-02 

I.  0023 6-31 


5.093  76  +  19 
5.93256+19 
5.91  1 36+19 
5.91956+19 
5.92516+19 
5.93636+19 
5.94636+ 1 9 
5.95216+19 
5 . 9  5  9  B  6  +  1 9 
5.96  746+  19 
5.97466+19 
5  .952  16  +  1.9 
5.95926+19 
5.99626+19 
6.03316+19 
5.03996+19 
6.01656+19 
6.02316+19 
6.02956+19 
6.03586+19 
6.04236+19 
6 • 04  8  76  +  19 
6.05426+19 


7.16936+19 
7.01886+19 
6.87356+19 
6.73346+19 
6.5980C+L9 
6.46736+19 
6.34086+19 
6.21866+19 
6.100  36  +  19 
5.98586+19 
6.87496+ 19 
5.76756+19 
5.66336+ 19 
5.56236+1 9 
5.46446+19 
8.36936+19 
5.2770C+19 
5.18746+1 9 
5 .10046  ♦  19 
5.31576+19 
4.93356+19 
4.85356+19 
4.77576+19 


